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ABSTRACT 



After a brief introduction to the quantum no-cloning theorem and its link with 

the linearity and causality of quantum mechanics, the conccpt of quantum cloning 
machines is sketched, foUowing, whenever possible, the chronology of the main re- 
sulta. The important classes of quantum cloning machines are reviewed, in partic- 
ular state-independent and state-dependent cloning machines. The l-to-2 cloning 
problem is then studicd from a formal point of vicw, using the isomorphism between 
completcly positive maps and operators, which leads to the so-called doublc-BcU 
ansatz. This also yields an cfhcient numerical approach to quantum cloning, bascd 
on semidefinite programming methods. The dcrivation of the optimal N-to-M uni- 
versal cloning machine in d dimensions is then detailed, as well as the notion of 
asymmetric cloning machines. In the second part of this revicw, the optical imple- 
mentation of cloning machines is considcred. It is shown that the universal cloning 
of photons can bc achievcd by paramctric amplification of light or by symmetriza- 
tion via the Hong-Ou-Mandel effect. The various experimental dcmonstrations of 
quantum cloning machines are reviewed. The cloning of orthogonally polarizcd 
photons is also considered, as well as the asymmetric and phase-covariant cloning 
of photons. Finally, the extension of quantum cloning to continuous variables is 
analyzed. The optimal cloning of coherent states of light by phase-insensitive ampli- 
fication is explained, as well as the experimental realization of continuous-variable 
quantum cloning with linear òptics, measurement, and feed-forward operations. 
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§ 1. Introduction and history 
1.1. THE NO-CLONING THEOREM 

The history of quantum cloning can bc traccd back to the controversial story of a 
paper by Herbert [1982] entitled "FLASH - A superluminal communicator based 
upon a new kind of measurement" . In this paper, submitted early 1981 to Founda- 
tions of Physics, Herbert was discussing an idealized làser gain tube which would 
produce, via sümulated emission, macroscopically distinguishable states of hght for 
an incoming singlo photon in any polarization state. The claim was that the noise 
in this process would, in principle, not prevent perfectly identifying the polarization 
state of the photon. This process would supposedly open the way to faster-than- 
light communication, a possibility with which any physicist feels uncomfortable 
since it violates causality. 

Today, morè than twenty years later, it is publicly known that GianCarlo Ghi- 
rardi and Asher Peres were requested to review this paper. The first of them 
recommended its rejection, based on the argument that the linear nature of quan- 
tum mechanics must prevent such a process to exist, see van der Merwe [2002]. The 
second referee wrote, see Peres [2002], that he had realized the paper was wrong, 
but neverthelcss recommended its publication bccausc he oxpccted that finding the 
error would raisc a considerable interest! Herbert 's paper was then published, and, 
funnily cnough, the prediction of Peres happened to be true. 

Soon afterwards, Woottcrs and Zurck [1982] published a paper in Nature, enti- 
tled "A single quantum cannot be cloned", which arrived at essentially the same 
conclusions as those drawn by Ghirardi in his anonymous referee report dated 
of April 1981, which itsolf was turncd into a paper two years later, see Ghirardi 
[1983]. Wootters and Zurek realized that, if one can build a "cloning machine" 
that produces several clones of the horizontal- and vertical-polarization states of an 
incoming photon, then circularly-polarizcd states cannot yield circularly-polarized 
clones. Instead, due to the linearity of quantum mechanics, one gets a linear super- 
position of vertically-polarized clones and horizontally-polarized clones. Indeed, if 
the cloning machine is such that 

\H) \C) - \H,H) \Ch), \V) \C) - \V,V) \Cv), (1.1) 

where \H) and |V) are horizontal- and vertical-polarization states of the original 
photon, \C) is the initial state of the cloning machine, and \Ch) and |Cy) are the 
(arbitrary) final states of the cloning machine, then the left and right circularly- 
polarized states, \L) = 2-'^/'^{\H) + i\V)) and \R) = 2-^l'^{\H) - i\V)), are trans- 
formed as 

|L) |C) ^ 2-i/2(|F, H) \Ch) + i\V, V) \Cv)) ^ \L, L) \Cl), 

\R) \c) ^ 2-i/2(|7ï, H) \Ch) - i\V, V) \Cv)) + \R. R) \Gr), (1.2) 

with |Cl) and \Cr) being (arbitrary) final states of the cloning machine. As a 
consequence, the cloning of circularly-polarized states fails, even in the special case 
where \Ch) = \Gv). This is the simplest explanation of what is known today as 
the quantum no-cloning theorem. 

Independently of this story, a related paper was published almost simultaneously 
by Dieks [1982] in Physics Letters, also showing that the "FLASH" proposal by 
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Herbert was flawcd. Hcrc, thc proof relies on the existencc of EPR statcs, see Ein- 
stein, Podolsky and Rosen [1935], which give rise to quantum correlations between 
spatially separated systems. If two photons are prepared in the EPR state 

|EPR) = 2-^/^(|íí, V) - I V, H)) (1.3) 

it is well known that mcasuring the hnear polarization of one of them in the 
horizontal- vertical basis allows one to immcdiately predict the outcome of a mca- 
surement of the linear polarization of the second one in the same basis, even if 
the measurement events are separated by a spacelike interval. For example, if the 
first photon is found to be in the \H) state, thcn the second photon will necessar- 
ily be observed in the \V) state. This property holds for any measurement basis. 
It had been realized since the early times of quantum mechanics that this prop- 
erty, callcd quantum entanglement, does not permit superluminal communication. 
Indeed, thc statistics of any measurement performed on one of the twin photons 
remains unchanged irrespectively of the measurement (or, more generally, the op- 
eration) applied on the second one. Dieks noticed, however, that if it was possible 
to clone perfectly one of the twin photons when the other one has been measured, 
then superluminal communication would become possible; hence, cloning must be 
impossible. 

Assume that Alicc measures the first photon either in thc horizontal-vertical lin- 
ear polarization basis or in the left-right circular polarization basis depcnding on 
whether shc wants to transmit a or a 1 to Bob. In the former casc, the second 
photon will be found by Bob to be in a balanccd mixturo of the \H) and \V) states, 
while in the second case it will be in a balanccd mixturc of the \L) and \R) states. 
These two mixtures are indistinguishable (they are characterized by the same den- 
sity operator, proportional to the identity J), which is why quantum mechanics is 
said to "coexist peaccfully" with spccial rclativity. However, if the second photon 
could be cloned perfectly, Bob would then get either a balanced mixturc of \H, H) 
and |V,y), or a balanced mixturc of \L^L) and \R,R). These mixtures being dis- 
tinguishable, Bob would have a way to infer Alicc's bit instantancously (with some 
error, but which can be made arbitrarily small as the number of clones increases). 
Dieks concluded from this paradox that such a cloning transformation cannot be 
consistent with quantum mechanics. 

It appears that the quantum no-cloning theorem is thus one of these scientific 
results that has been rediscovered several times, at least by Dieks, Ghirardi, Woot- 
ters, and Zurek. Actually, it can be argued that it was alrcady implicitly used 
by Stephen Wiesner in his famous paper entitled "Conjugate coding" written in 
the 70's but only published in Wiesner [1983], which is sometimes considered as 
the founding paper of quantum Information theory. In some sense, thc no-cloning 
theorem was already intrinsically contained in the roots of quantum mechanics and 
is thus trivial; but, on the other hand, its discovery has contributed to revisiting 
quantum mechanics in an Information thcorctic language, which has had a decisive 
influence on the dramàtic development of quantum Information science over the 
past twenty years. 
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1.2. BEYOND THE NO-CLONING THEOREM 

Soon aftcr thc publication of thc quantum no-cloning thcorcm, anothcr paper ap- 
peared in Nature, written by Mandel [1983]. In this paper, entitled "Is a photon 
amplifier always polarization dependent?", Mandel drew attention to the physi- 
cal origin of thc impossibility of making a pcrfcct aniplifying apparatus for light, 
namely spontaneous emission. He showed that, if the amphfier is a single 2-level 
atom with a dipole moment /í, then the amplification of an incoming photon of 
polarization vector e dcpcnds on thc scalar product bctwecn /i and e. If thc polar- 
ization vector e of the incoming photon is parallel to the dipole moment /i, then 
the state |l)e will, after some interaction time, evolve into a state containing the 
dcsircd two-photon state |2)t duc to stimulated emission. On the contrary, if e is 
orthogonal to n, then the two-photon component of the resulting state corresponds 
to |l)e|l)e, where e is a polarization vector orthogonal to e. This is due to sponta- 
neous emission, which spoils the amplification since one of thc two photons has thc 
wrong polarization e. In other words, with such a simple one-atom amplifier, the 
final state depends on the polarization of the incoming photon. 

Interestingly, Mandel noticed that if we consider a more elaborate amplifier made 
of two such àtoms with orthogonal dipole moments (/Xi and 112), it becomes never- 
theless possible to amplify the photon independently of its polarization, although 
this process suífers from thc unavoidable noise originating from spontaneous emis- 
sion. Assuming that the two àtoms interact similarly with the incoming photon, 
one understands intuitively that if one atom amplifies the photon "well" (when e 
is close to /ii), then the second atom amplifies it "poorly" (because e is then ap- 
proximately orthogonal to 112)- The balance between these two efïects results in an 
amplification that does not depend on e. By filtering out the resulting two-photon 
component, one gets 

|1), |0), - l |2),(2| ® |0),(0| + i |1),(1| |1),(1| (1.4) 

irrespective of e. In some sense, the perfcct cloning of polarization via stimulated 
emission works with probability 2/3, while spontaneous emission blurs the polar- 
ization with probability 1/3. 

Mandel's paper remained mostly unnoticed and, remarkably, one had to wait 
more than ten ycars before the notion of quantum cloning machine, which was 
implicitly contained in this paper, became popular. In a seminal paper, Buzck and 
Hillery [1996] realized that, although perfect quantum cloning is ruled out by thc 
no-cloning principle, some imperfect cloning may be possible. They found out that a 
qubit (2-level quantum system) that is in an unknown state can be approximately 
duplicated, resulting in two pretty good clones of the original state. This result 
holds in full gcnerality, rcgardless of the physical variable carrying the qubit, so 
it goes much beyond the polarization-independent amplification of a single photon 
considered bc^fore. This paper had a considerable impact at that time because 
quantum Information thcory was born, and it had bcen realized how fruitful it is 
to investigate quantum mechanics using an Information language. 

Consider a qubit in the state = a\0) + f3\l), where |0) and |1) form an or- 
thonormal basis of the Hilbert space, while a and /3 are arbitrary complex numbers 
satisfying jaj^ + |/3|^ = 1. Buzek and Hillery [1996] addressed the foUowing formal 
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Figure 1: Quantum cloning machinc. The input of thc transformation consists of 
the state \tp) to be cloned and the initial state of the cloning machine |C). The two 
clones are contained in outputs A and B, while C refers to an ancilla. 
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problem: find a transformation acting on an original qubit in state \ip) together 
with an auxiliary system (commonly viewed as the cloning machine itself) that 
produces two clones with the same fidelity and is state-independent, or universal, 
li the cloning machine is initially put in state |C), then 

|0)|C)-|Eo), (1.5) 

with the final states |Eo) and |Ei) belonging to the product Hilbert space TÍa ® 
Hb^ "Hc, where Ha and Hb denote the spaces of the two clones (called A and B) 
and He denotes the space of thc cloning machine C, see Fig. 1. By linearity, an 
arbitrary qubit state \tp) is cloned as 

IV) |C) ^ a|So)+/3|Ei) = |E). (1.6) 

The fidelity of the clones, which measures the overlap between the input state and 
each clone, is given by 

fAW = (VITrsc(S)IV) , /b(V) = ii^lTrAci^m , (1.7) 

where Tr denotes the trace and S = |E)(E| is a short-hand notation for the density 

operator of a pure state. Buzek and Hillery [1996] showed that, under the constraint 
that /a(V') = /b(V') is independent oíijj, quantum mechanics permits the existence 
of a cloning transformation which achieves a fidelity as high as 

= 5/6 ~ 0.833 (1.8) 

This transformation, which is called a quantum cloning machine, is given by 

|0)|C) ^ |Eo) = 7273 |00)ab |0)c + 073 |l)c, 

|1)|C7) ^ |Ei) = y273|ll)AB |l)c + yï73|*+)As |0)c, (1.9) 

where |*+) = 2-V2(|oi) + 110)) is one of the Bell states, while |0)c and |l)c denote 
two orthogonal states of the cloning machine. It is easy to check, by tracing over 
the cloning machine, that the two clones of an input state |0) are left in the joint 
state 

PAB = Trc(E) = ^ |00)(00| + ^|*+)(*+| (l·lO) 
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which is equivalent to Eq. (1.4) given the bosonic statistics of photons. More gen- 
erally, if the input state is jí/;), the first term of the right-hand side of Eq. (1.10) 
becomes a projector onto while the second term is some (í/;-depending) 

maximally-entangled statc. Therefore, by tracing over one of the clones, the result- 
ing state of the other clone is 



with I denoting the identity opcrator, confirming that the two clones are left in 
the same state. Tlicy can bc vicwcd cach as cmcrging from a quantum depolarízing 
channel: they are found in the right state \tp) with probability 2/3, while they are 
replaced by a random qubit 1/2 with probability 1/3. 

Soon aftcr tlic publication of tliis paper, it was provcd that this niachinc is 
actually the optimal universal cloning machine, that is, the highest fidelity of cloning 
permitted by quantum mechanics is indeed 5/6, see Bruss, DiVincenzo, Ekert, 
Fuchs, Macchiavello and Smolin [1998]. 

This discovery by Buzek and Hillery [1996] triggered an immense interest and 
initiated an entire subfield of quantum Information science devoted to quantum 
cloning. In particular, further studies addresscd cloning in dimensions largcr than 
2, state-dependent cloning (considering a restricted set of input states), the so- 
called N-to-M cloning (where one produces M identical clones out of N identical 
replicas of the original), the asymmetric cloning (where the clones havc uncqual 
fidelities), the cloning of orthogonal qubit states, the cloning of continuous- variable 
states (such as coherent states), the economical cloning (where no ancillary space 
is necessary), the probabilistic cloning (which is not dcterministic, i.e., it does 
not succeed with probability 100%), or even the cloning of quantum entanglement 
(instead of quantum states). These numerous results will be reviewed in Section 2. 

Aside from its utmost importance for the foundations of quantum mechanics, the 
study of quantum cloning has drawn a lot of interest probably also because it is 
closely connected to quantum key distribution (QKD). Indeed, in many cases, the 
cloning machine is known to bc the most powerful cavesdropping stratcgy against 
QKD protocols: the eavesdropper duplicates the quantum state and sends one clone 
to the authorized party, while keeping the second clone for later measurement. The 
characterization of cloning machincs is therefore crucial to assess the sccurity of 
these QKD protocols (this particular connection is out of the scope of the present 
review paper, and will not be discussed any further). 

1.3. QUANTUM CLONING WITHOUT SIGNALING 

Before entering the detailed study of quantum cloning, it is interesting to backtrack 
for a moment and further discuss the proof of the quantum no-cloning theorem 
based on a pair of entangled photons due to Dieks [1982]. As explained earlier, if 
Alice measures her photon either in the horizontal- vertical linear polarization basis 
or in the left-right circular polarization basis, and if Bob is able to clone his photon 
perfectly, then he obtains two distinguishable 2-photon mixtures, which apparently 
makes superluminal signaling possible. 

A natural idea, due to Gisin [1998], is to assume that Bob's cloning machine must 
necessarily introduce some intrinsic noise, and determine the minimum amount of 



p^ = TrBc(S) = ^ + 



PB=TV^c(S) = ^lV)(^| + ^/, (1.11) 
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noise that must be added so that causality ceases to be violated. Rcmarkably, it so 
happens that the minimum noise needed to comply with causality exactly coincides 
with that of the optimal universal cloning machine. In other words, the upper 
bound on quantum cloning can be derived from simple principies. 

As shown earlier, the two clones of the universal machine emerge cach from a 
quantum depolarizing channel, see Buzek and Hillery [1996]. This channcl can be 
interpretcd as giving rise to a shrinking of the vector representing the qubit state 
in the Bloch sphere. Using the Bloch rcpresentation 

p = (/ + m • a)/2, (1.12) 

where m is a vector isomorphic to state p and cr = {ax,cry, (Jz) is the vector of Pauli 
matrices, we see from Eq. (1.11) that a state associated with m yields two clones 
which are in a state associated to 2m/3, independently of the orientation of m. 
Therefore, this universal cloning machine is sometimes also termed as isotropic. 

Following Gisin [1998], consider that the (pure) state of the original qubit is 
associated with the (unit-norm) vector m, and Ict us rcstric;t our search to cloning 
machines that are symmetric and isotropic, that is, the clones are in the states 

PA(m) =pB(m) = (l + r/m·a)/2, (1.13) 

where ry is an unknown "shrinking factor" (0 < r/ < 1). It is easy to check that r] is 
related to the fidelity by 

/^(m)=/B(m) = (l + r?)/2. (1.14) 

Using Eq. (1.13), the Hilbert-Schmidt decomposition of the joint state of the two 
clones can be written as 

p^s(m) = {l + • (7 7 + / 7?m • cr + ^íj.fcCTj ® a·fcj/4, (1-15) 

j-.k 

where the matrix tj^k measures the quantum correlations between the clones. Gisin 
[1998] went on deriving constraints on tj^k that result from covariance and causality. 
The covariance property (which will be explaincd in detalls latcr on) mcans, phys- 
ically, that rotating the original qubit around, say, the z-axis before cloning must 
be equivalent to cloning the original qubit and then rotating each of the two clones 
by the samc amount around the z-axis. Following Dieks' argument, the causality 
condition is taken into account by imposing that 

PAB(mi) + pAB(-mi) = pAB(m2) + pAB(-m2) (1.16) 

which expresses the fact that the two-clone states corrcsponding to two indistin- 
guishable mixtures of input states, {mi,— mi} and {m2, — m2}, are themselves 
indistinguishable. Putting all these conditions on tj^k together, one can show 
that the maximum value of 77 that preserves the positivity of the two-clone state, 
Pab{t^) > 0, is ry = 2/3; hence fA = fs = 5/6. This provides an alternate proof of 
the optimality of the qubit universal cloner of Buzek and Hillery [1996]. 

To be complete, Ict us mention that such a use of the no-signaling condition 
has been criticized in Bruss, D'Ariano, Macchiavello and Sacchi [2000], the ar- 
gument being that the linearity and trace-preservation properties of the cloning 
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map (which, combined, imply the no-signaling condition) are not sufíicicnt, strictly 
speaking, and need to be supplemented with the complete positivity condition in 
order to bound the cloning fidelity. This simple technique, howevcr, has proved 
to be successful to recover conditions on probabiHstic cloning, see Hardy and Song 
[1999], on asymmetric universal cloning, see Ghosh, Kar and Roy [1999], or even to 
íind a new class of real cloning machines, see Navez and Cerf [2003]. 

1.4. CONTENT OF THIS REVIEW PAPER 

The rest of this review paper will be devoted to the study of quantum cloning 

machines, as well as their optical realization. Let us sketch the content of the 
next Sections. Section 2 will provide an overview of the chronology of the main 
papers that have been written in this context, focusing on the results but skip- 
ping the dcrivations. The numerous classes of quantum cloning machines will be 
presented (universal cloners, Pauli or Heisenberg cloners, phase-covariant cloners, 
Fourier-covariant cloners, group-covariant cloners, real cloners, entanglement clon- 
ers, continuous- variable cloners, probabilistic cloners, or cconomical cloners). 

In Section 3, we will consider the issue of quantum cloning from a formal point of 
view, based on the description of the associated completely-positive (CP) map and 
the notion of covariancc. This study will be rcstrictcd to l-to-2 cloning, and will 
focus on the isomorphism between CP maps and operators. It will be shown that 
finding the optimal cloning map reduces to a semidefinite programming problem, 
which can be cfficicntly solvcd by numcrical mcthods. It will also be shown that 
the unitary realization of a cloning map based on the "double-Bell' ansatz provides 
a simple and efRcient tool to investigate cloning analytically. Some examples of 
d-dimcnsional l-to-2 cloners will be providcd. 

In Section 4, this formal study will be extended to N-to-M cloning machines in 
d dimensions, but will be restricted the case of universal cloning. The derivation 
of the optimal cloning transformation as well as the optimality proof will be dc- 
tailed. In addition, the extension to asymmetric cloning machines and the notion 
of universal-NOT gate will be discussed. The reader who is mainly interested in 
the optical realization of cloning machines and not so much in their theoretical 
derivation may skip Sections 3 and 4, and immediately go to the next Sections. 

In Section 5, the optical implementation of the universal quantum cloning ma- 
chines will be analyzcd in detalls. The cloning experiments relying on stimulated 
parametric downconversion will be described first, foUowed by those relying on 
the symmetrization that can be obtained with a Hong-Ou-Mandel interferometer. 
Next, the optical realization of (universal) asymmetric cloning machines will be 
discussed, as well as the (universal) cloning of a pair of orthogonal qubits. 

In Section 6, the phase-covariant cloning machines will be developed for qubits 
as well as d-dimensional systems, in a l-to-2 or iV-to-Aí configuration. The experi- 
mental realization of phase-covariant cloning for photonic qubits will be described. 

In Section 7, the generalization of quantum cloning to states bclonging to an 
infinite-dimensional Hilbert space will be considered. In particular, the cloning of 
coherent states of light by phase-sensitive amplification will be explained, as well 
as the experimental realization of continuous-variable cloning using linear òptics, 
measurement, and feed-forward. The cloning of a finite-width distribution of co- 
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herent states will bc analyzcd, as wcll as thc cloning of a pair of conjugate coherent 
States. Finally, the conclusions will be drawn in Section 8. 

§ 2. Overview of quantum cloning machines 
2.1. UNIVERSAL CLONING MACHINES 

This Section will be devoted to a summary of the various cloning machines that 

have been introduced in thc litcrature, foUowing the chronology as well as possible. 
Soon after the universal quantum cloning machine was discovered by Buzek and 
Hillery [1996], the question arose whether this machine was optimal. As already 
mentioned, this cloning machine is required to be symmetric, that is, thc 2 clones 
must have equal fidelities /a(V') = ÍbÍÍ^), V-í/;. In addition, it must bc universal (or 
state- independent), which means that all states are cloned with thc samc fidelity, 
independent of ip. It was proven by Bruss, DiVinccnzo, Ekcrt, Fuchs, MacchiavcUo 
and Smolin [1998] that it is indeed the optimal symmctric universal duplicator for 
qubits, so that / = 5/6 is indeed the highest fidelity allowcd by quantum mechanics 
in this casc. In the same paper, the concept of optimal state- dependent cloning 
machines was also introduced, that is, transformations that optimally duplicate 
orily a particular subsct of the input states. Almost simultaneously, Gisin and 
Massar [1997] introduced thc concept of N-lo-M quantum cloning machines, which 
transform N idcntical rcplicas of an arbitrary state, jV')®'^, into M > N identical 
clones. They were able to prové for low N that the optimal universal N-to-M 
cloning of qubits is characterized by the fidelity 

M{N + 1) + N 
M{N + 2) 

Incidentally, this confirms the optimality of the l-to-2 universal cloning machine 

with fidelity = "V^- Thc quantum nctwork that rcalizcs this l-to-2 imivcrsal 

cloning of qubits was described by Buzek, Braunstein, Hillery and Bruss [1997], 
while it was extended to 1-to-M universal cloning in Buzek and Hillery [1998b]. 
Note also that when thc number of clones M increases for fixcd iV, thc cloning 
fidelity decreases. This can simply be interpreted as a spreading of quantum Infor- 
mation over more clones. In the limit M — > oo, the cloning transformation tends 
to a mcasurcmcnt, which confirms that thc optimal (state- independent) estimation 
of the state of N identical qubits has a fidelity 

^univ ' /Q 9\ 

/JV^OO — N + 2 

as originally dcrivcd in Massar and Popcscu [1995]. 

Then, in early 1998, the extension of quantum cloning machines to higher- 
dimensional spaces was considered independently by Buzek and Hillery [1998a], 
Ccrf [1998], Wcrncr [1998]. Thc form of thc optimal universal l-to-2 cloncr in di- 
mension d was conjectured by Buzek and Hillery [1998a], Cerf [1998], while the 
derivation and full optimality proof of the universal d-dimensional N-to-M cloner 
was givcn by Wcrncr [1998], Kcyl and Wcrncr [1999]. Thc optimal fidelity of the 
universal l-to-2 cloner of d-dimensional states (or qudits) was shown to be 



j-univ _ -l; -r i» , , 
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while, for arbitrary N and M > N,\i is 

.univ ... M(iV+i) + (rf-i)Ar 

In Cerf [1998], the cloning of d-dimensional systems was actually investigated 
in a more general setting: a large class of symmetric or asymmetric, universal or 
statc-dcpcndcnt, l-to-2 cloning machincs was introduccd in arbitrary dimcnsion d. 
The optimality of this class of cloners was only conjectured, but, in the special 
case of a symmetric and universal cloner, Eq. (2.3) was also derived. For the set of 
asymmetric imiversal l-to-2 cloning machines, the balance between the fidelity of 
the two cloncs 

fT'^id) =T1A + {1- VA)/d, =71B + {1- VB)/d, (2.5) 

was characterized by the simple relations 

VA = l-a\ VB = l-í3\ a' + ^+(3' = l, (2.6) 

where tja and tjb are the "shrinking" factors associated with the clones {rj is the 
probability that the input state emerges unchanged at the output of the quantum 
depolarizing channel). Here, a and /? are positive real variables. It is instructive 
to notice that in the limit d ^ oo, the cloning of quantum Information resembles 
the distribution of a resource that can strictly not be shared: the probability that 
IV') is found in one clone is complementary to the probability that it is found in the 
second clone, that is, ija + rjB = I- 

The quantum circuit for the asymmetric universal cloning of qubits was displayed 
in Buzek, Hillery and Knight [1998], while it was generalized to d-dimensional 
symmetric or asymmetric cloning in Braunstein, Buzek and Hillery [2001]. 

Finally, even more general quantum cloning machines were obtained in the special 
case of qubits (d = 2) in an independent work by Niu and Grifíiths [1998]. There, 
the l-to-2 asymmetric and state-dependent cloning of a qubit was investigated in 
full generality, and, in particular, formulas (2.6) were recovered for d = 2 without 
any assumption. 



2.2. PAULI AND HEISENBERG CLONING MACHINES 

The results of Cerf [1998] were later expanded in Cerf [2000a] for the case of qubits, 

and Cerf [2000b] for the casc of d-dimcnsional systems. The specificity of the ap- 
proach to quantum cloning underlying these papers is that one considers the cloning 
of a system that is initially maximally entangled with another system instead of 
the cloning of a purc state. This second system acts as a "rcference" by kccping a 
memory of the original state after the cloning has been achieved. The final state of 
the "reference" , the two clones, and the cloning machine then fully characterizes the 
cloning transformation, as a consequence of the isomorphism between complctely 
positive (CP) maps and operators (this will be explained in detalls in Section 3.1). 
By choosing an appropriate form for this final state, one generates a large class of 
quantum cloning machines. 

For qubits {d = 2), this class corresponds to the so-called Pauli cloning machines, 
whose clones emerge from two - possibly distinct - Pauli channels. In a Pauli 
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channel, the input qubit undergoes one of the thrcc Paiili rotations {(Jxi'^yi'^z} 
or the identity / with respective probabilities {PxTPy,Pz, 1 -~ Px ~ Py ~ Pz}- For 
example, it was shown that the whole class of symmetric PauH cloning machines 
corresponds to PauH channels with probabilities Px = x"^, Py = U^, and Pz = z'^·, 
with x, y, z satisfying the condition 



The action of these PauH cloners is easy to understand knowing that, if the original 
qubit is in an eigenstate of ax, namely (|0) ± |l))/·\/2, then it is rotated by an 
angle tt around the y-axis (z-axis) under ay (az) while it is left unchanged (up to 
a sign) by ct,,.. Thcrefore, the cloning fidelity of the eigenstates of ax 1 — py — pz- 
Similarly, the eigenstates of ay, namely (|0) are cloned with the fidelity 

1—Px—Pz, while the eigenstates of a^, namely |0) and |1), are cloned with the fidelity 
1 — Px — Py The universal l-to-2 symmetric cloning machinc simply corresponds 
to Px = Py = Pz = 1/12. Note that these Pauli cloning machines appear to be a 
special case of the state-dependent cloning transformations considered in Niu and 
Grifíiths [1998]. 

These considerations can be extended to d dimensions in order to obtain the set 
of so-called Heisenberg cloning machines, whose clones emerge from two - possibly 
distinct - Heisenberg channels. In a Heisenberg channel, the d-dimcnsional input 
state undergoes, according to some probability distribution, one of the cP error 
operators í^m.n (with < m,n < d — 1) that form the discrete Weyl-Heisenberg 
group. It can bc shown that the probability distribution of the errors for 

the first clone is dual, under a Fourier transform, to that of the second clone. 
This corroborates the fact that if one clone is close-to-perfect (its associated error 
distribution is peaked), then the second clone is very noisy (its associated error 
distribution is flat). More precisely, this fidelity balance between the two clones 
can be shown to result from a no-cloning uncertainty principle, akin to Heisenberg 
principio, sce also Cerf [1999]. The quantum circuit realizing these Heisenberg 
cloning machines was displayed in Braunstein, Buzek and Hillery [2001]. 

Recently, the optimality of this entire class of (Pauli or Heisenberg) quantum 
cloning machines has been rigorously proven by Chiribella, D'Ariano, Perinotti 
and Cerf [2005] in the foUowing sense: under some general invariance conditions, 
the cloners of this class coincide with all the extremal cloners. Therefore, for a given 
(invariant) figure of mèrit, it is sufficicnt to search the optimal cloner within this 
class to be guaranteed that the found solution is the global optimal cloner. 

2.3. PHASE- AND FOURIER-COVARIANT CLONING MACHINES 

In 2000, an important class of statc-dcpcndcnt qubit cloning machines. named 
phase-covariant cloning machines, was introduced by Bruss, Cinchetti, D'Ariano 
and Macchiavello [2000] . It is defined as a transformation that clones all the balanced 
superpositions of basis states with the same (and highest) fidelity. These states 



x'^ + y'^ + z^ + xy + xz + yz 



1/2 



(2.7) 



|V') = (|0)+e^1l))/^/2 



(2.8) 



with (j) being an arbitrary phase, are located on the equator of the Bloch sphere. 
The optimal cloner also fulfills the covariance condition with respect to the rotation 
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of (f>, that is, cloning the rotatcd original qubit is equivalent to cloning the original 
qubit foUowed by a rotation of each of the clones. The optimal phase-covariant 
symmetric l-to-2 cloner was found to have a fidelity 



which is higher than that of the corresponding universal cloner, fi^2('^) = 5/6. In 
contrast, the resulting fidelity for the states |0) and corresponding to the poles 
of the Bloch sphere, is equal to 3/4, which is lower than /"i^'JÍ^)- In some sense, it 
is possible to bettcr clone some restricted set of states (the equator) at the expense 
of a worse cloning of sonic other states (near the poles) . 

Interestingly, this phase-covariant cloner can be viewed simply as a special case 
of the Pauli cloncrs, sce Ccrf, Durt and Gisin [2002]. If wc take x = y = l/\/8 
and z = 1/2 - l/\/8, which satisfies Eq. (2.7), we indeed recover the same cloner: 
the eigenstates of ax are cloned with fidelity í — Py — Pz = 1/2 + l/'v/S, while 
the eigenstates of CTj, are cloned with the fidelity l - Px - Pz = 1/2 + l/VS. It 
was observed in Cerf, Durt and Gisin [2002] that imposing these 4 states lying 
symmetrically on the equator to be cloned with the same fidelity results in the 
phase-covariant cloner, which actually gives the same fidelity for all states on the 
equator (the deep reason for this equivalence was found in Chiribella, D'Ariano, 
Perinotti and Cerf [2005]). Finally, we verify that this Pauli cloner clones the 
eigenstates of cr^ with a lower fidelity 1 — Px — Py = 3/4. 

One can summarize the results on qubit cloning machines by noting that the 
eigenstates of the three Pauli matrices play the role of three mutually unbiased 
(MU) bases for qubits (MU bases are such that the modulus of the scalar product 
of any two states taken from distinct bases is \/\/d, with d being the dimension). 
One can thus define three genèric classes of qubit cloning machines, namcly, the 
universal cloner (which can be obtained by imposing the states of 3 MU bases to be 
cloned with the same and highest fidelity), the phase-covariant cloner (if the states 
of only 2 MU bases are cloned equally), and some particular Pauli cloner (if the 
states of all 3 MU bases are cloned with unequal fidelities). The cloning of qubits 
having been essentially covered, it became natural to turn to the state-dependent 
cloning of qutrits {d = 3). 

In Cerf, Durt and Gisin [2002] , four kinds of Heisenberg cloning machines were 
defined for qutrits, depending on whether four, three, two or none of the MU bases 
are cloned with the same fidelity. If none of the MU bases are cloned with equal 
fidelities, one has a particular Heisenberg cloning machine. On the contrary, if all 
four MU bases are cloned with the same fidelity, one recovers Eq. (2.3) for d = 3 in 
the case of symmetric cloning, that is, the qutrit universal cloner with fidelity 



If three MU bases are requested to be cloned with the same fidelity, one gets the 
so-called double-phase covariant qutrit cloner, with a fidelity 



(2.9) 



/nj(3) = 3/4. 



(2.10) 



7^2(3) = 



5 + VI7 
12 



~ 0.760 



(2.11) 



slightly higher than fi^2{'^)- This cloner, which was independently derived in 
D'Ariano and Lo Presti [2001], has the property that it clones with the same and 
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highest fidelity all the balanced superpositions 

|V) = ^(|0)+e*|l)+e^^^|2)) (2.12) 

for arbitrary phases 4>i and 02 (it is also covariant with respect to both (pi- and (/>2- 
rotations). This can be understood by noting that if we complete the computational 
basis with any triplct of bascs in ordcr to makc 4 MU bascs, thcsc 3 bases only 
consist of balanced superposition states. In analogy with the qubit case, it then 
appears that imposing these 3 bases to be cloned with the same (and highest) 
fidelity rcsnlts in a cloning machine that clones all states (2.12) equally well, that 
is, the double-phase-covariant cloner. 

Finally, we may impose that two MU bases that are dual under a Fourier trans- 
form arc cloned with the sanic (and highest) fidelity, the othcr two bcing also 
cloned with an equal (albeit lower) fidelity. For cxample, the computational ba- 
sis {|0), |1), |2)} and the dual basis = 3-^/2 ^^^^ i^'^lfc) with j = 0,1,2 and 
-y = form such a pair of MU bascs. Wc then gct the so-callcd Fourier- 

covariant cloner for qutrits, see Cerf, Durt and Gisin [2002], with a fidelity 

_^Fourier(3) ^ ^ ^ ^ 0.789 (2.13) 

which is evcn liiglicr tlian ,fi!ï^2Í'^) cxpcctcd sincc the considcrcd set of input 
states is smaller than for the double-phase-covariant cloner. This cloner is covariant 
with respect to a Fourier transform, hence it clones two Fourier-conjugate bases with 
the same fidelity. 

Note that, except in dimension 2, any two MU bascs cannot always be mapped 
onto any other two MU bases, so that the Fouricr-covariant cloner is not the unique 
transformation that clones equally well two MU bases. Indeed, in Durt and Nagler 
[2003], it was shown that, in dimension 4, the cloner for two MU bases conjugate 
under a Fourier transform differs from the cloner for two MU bascs conjiigate under 
a double Hadamard transform. In the spccial case of qubits {d — 2), however, all 
pairs of MU bases are unitarily equivalent, so that the Fourier-covariant and phase- 
covariant cloners coincide, fl^2^'^^{2) = /f^2(2)· 



2.4. GROUP-COVARIANT CLONING MACHINES 

In D'Ariano and Lo Presti [2001], a general method for optimizing the group- 

covariant cloners was dcrivcd. Morc specifically, thcy considcrcd the optimal 
cloning transformations that are covariant under a proper subgroup Í2 of the uni- 
versal unitary group U (d). For example, the universal qubit cloner is covariant with 
respect to U{2), whilc the phase-covariant qubit cloner is covariant with respect to 
U{1). They used this technique to derive the symmetric double-phase-covariant 
cloner for qutrits corresponding to Eq. (2.11), as well as the l-to-3 symmetric 
phase-covariant cloner for qubits, associated with the fidelity 

/f^3(2) = ^ ~ 0.833. (2.14) 

Owing to the complexity of the group-theoretical parametrization of CP maps un- 
derlying this technique, its applicability seems rather limited. Nevertheless, using 
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another method, Fan, Matsumoto, Wang and Wadati [2001] were able to derive the 
optimal l-to-AÍ symmetric phase-covariant cloning of qubits, yielding the fidelity 



1 JM(M+2) 

fvo (2.) 2 + Meven, 

More reccntly, D'Ariano and Macchiavcllo [2003] succeedcd in applying this theory 
of group-covariant cloning in ordcr to confirm Eq. (2.15), as wcll as to find a gen- 
eral expression for /^"l^^^(2) and tlic associated N-to-M cloner. This cxpression, 
which was partly conjccturcd in Fan, Matsumoto, Wang and Wadati [2001], is quite 
complex, and depcnds on whcthcr N and M havc the samc parity. It was noticed 
that, if the parities do not match, thcn the cloner that optimizes the fidelity of each 
of the clones does not coincide with the optimal cloner with respect to the global 
fidelity (measuring how well the joint state of the clones approximates {ip)'^'^'^ , if 
is the state of the original). In the case of qutrits {d = 3), D'Ariano and Macchi- 
avcllo [2003] also found the optimal 1-to-M symmetric phase-covariant cloner. The 
expression for its fidelity /f^jv^(3) is rather complex, and depends on M modulo 3. 

2.5. HIGH-D STATE-DEPENDENT CLONING MACHINES 

In parallel with this series of results on group-covariant cloning involving sevcral 
originals and clones but in low dimensions, both the phase-covariant and Fourier- 
covariant l-to-2 cloning machines were extended to arbitrary dimensions d. In 
Cerf, Bourennane, Karlsson and Gisin [2002], the d-dimcnsional symmetric Fourier- 
covariant cloner was derived, and shown to be characterized by the fidelity 

/fofer(^) = i + -L. (2.16) 

It clones equally well two MU bases that are conjugate under a Fourier trans- 
form, such as the computational basis {|0), • • • [d — 1)} and the dual basis — 
^-1/2 J2'l^Q 7^*^!^) with j = 0, • • • d — 1 and 7 = e^'^*/''. The asymmetric Fourier- 
covariant cloners were also characterized in the same paper. 

Thcn, in Fan, Imai, Matsumoto and Wang [2003], the rf-dimensional symmetric 
multi-phase-covariant cloner was derived, giving the fidelity 



fPc , d-2 + ^d^+M-4) 



A-2W = ^ + ^'-''^ 

It clones with the samc; (and highest) fidelity all balanced superpositions of the 
States of the computational basis, with arbitrary phases. This result was indepen- 
dently derived in Lamoureux and Cerf [2005], Rezakhani, Siadatnejad and Ghaderi 
[2005], where it was also extended to asymmetric cloners. Note also that the role of 
multi-phase-covariant cloners in the context of entanglement-bascd QKD protocols 
was first studied in Durt, Cerf, Gisin and Zukowski [2003] for qutrits, then in Durt, 
Kaszlikowski, Chen and Kwek [2004] for d-dimensional systems. 

2.6. CLONING A PAIR OF ORTHOGONAL QUBITS 

Another possible variant to the problem of cloning was studied by Fiuràsek, Iblisdir, 
Massar and Cerf [2002] , who introduced universal cloning machines that transform 
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2 qubits that arc in an antiparallcl joint statc j-í/i) j-í/'-'-) into M cloncs of with 
{ip\ip-^) = 0. It was proven that, when M is large enough, such a cloner outperforms 
the Standard 2-to-M cloner. One has the fidehty 



.univ _ 1 VWT2)7(3M) 

Ji.i^Ml^J - 2 2 ^ ^ ^ 

which is greater than /™5^(2) for M > 6. In some sense, it is better to replace one 
of the two original statí^s by its orthogonal state | ■)/'"'") if the goal is to produce 
M > 6 clones. This effect can be understood at the limit M oo, 

/iTZoo(2) = ^ + ^ ^ 0.789, (2.19) 

that is, for the optimal measurement of a pair of antiparallel qubits. Indeed, it had 
been noticcd carUcr by Gisin and Popescu [1999] that mcasuring 1-0) IV-"^) yiclds 
more Information than measuring with /^toÍ^) ~ '^/^· interpretation 

of this property lies in the dimension of the Hilbert space spanned by |V'"'")> which 
is 4, while 1^)®^ only spans the 3-diinensional symmetric subspace of 2 qubits, W+. 

2.7. ENTANGLEMENT CLONING MACHINES 

Another problem, related to quantum cloning, has been investigated in Lamoureux, 
Navez, Fiuràsek and Cerf [2004]. There, it was shown that the amount of entangle- 
ment contained in a two-qubit state cannot be cloned exactly, in analogy with the 
impossibility of cloning the state itself. If a cloning machine is devised that produces 
maximally-entangled clones for maximally-cntangicd rjubit pairs at the input, tlicn 
it cannot yield unentangled clones for all product states at the input. Neverthe- 
less, the approximate cloning of entanglement is very well possible. In Lamoureux, 
Navez. Fiuràsek and Cerf [2004], a class of l-to-2 entanglem,ent cloning machines 
was defined, which are universal over the set of maximally-entangled two-qubit 
states. The symmetric cloner of this class provides two clones of all maximally- 
entangled two-qubit states with the optimal fidelity 

/i«(2 X 2) = ^ 0.717 (2.20) 

corresponding to an entanglement of formation 0.285 c-bits. In contrast, all product 
states are transformed into unentangled clones. This was recently extended to the 
cloning of entanglement for d x d-dimensional systems in Karpov, Navez and Cerf 
[2005] . The fidelity of the optimal symmetric entanglement cloner that is universal 
over the set of maximally-entangled d x d-dimensional states is 



2.8. REAL CLONING MACHINES 




(2.21) 



Still another class of d-dimensional l-to-2 cloners was introduced in Navez and Cerf 
[2003] , named real cloning machines. It is defined as a transformation that clones 
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all real superpositions of tl·ie computational basis states with the same (and highest) 
fidelity. The optimal l-to-2 symmetric real cloner in dimension d was shown to have 
the fidelity 

, = i + 2-d + W-+4d + 20 
Ji^2\ ! 2 4(d + 2) ^ ' 

Note that in dimension d = 2, the set of real states forms a circlc in the Bloch 
sphere which is unitarily equivalent to the equator, so that we have fl^\{2) = 
/P^2(2) = /S2"'"'(2) = 1/2 + 1/V8. For any dimension d > 2, one has 

fiT2{d) < ffUd) < fi%{d) < ff^T'^id). (2.23) 

Note that, for d = 4, we have the identity /n'2(4) = f^'^2'^^{2 x 2). This comes 
from the fact that the set of maximally-entangled two-qubit states is isomorphic to 
the set of four-dimensional real states. 



2.9. HIGHLY ASYMMETRIC CLONING MACHINES 

In Iblisdir, Acín, Cerf, Filip, Fiuràsck and Cisin [2005] the conccpt of multipartitc 
A^-to-M cloning machines (with M > 2) was introduced. These machines are 
highly asymmetric in the sense that they produce M clones of unequal fidelities. 
A very general group theorctical approach to the construction of the multipartitc 
asymmetric cloning machines for qubits was then presented in Iblisdir, Acín and 
Cisin [2005]. It was applied to several particular examples such as the asymmetric 
1 ^ A'' + 1 cloning machine, which produces two kinds of clones, one clone with 
fidelity and N clones with fidelity f^. The optimal fidelities rc^ad 



fA 



= 1 - 



/f- 



1 



1 



.JV+l 



(2.24) 



where x G (0, 1) parametrizes the class of the optimal 1 — > A'+l asymmetric cloners. 
Note that the formula (2.24) holds only for A > 1. It also was conjecturcd, based 
on exact analytical calculations for low A^, that the optimal A^ ^ A^ + 1 asymmetric 
cloner which produces, from N replicas of a qubit, N clones with fidelity and a 
single clone with fidelity achieves 



'N+l 



f. 



N{N + 2)'"' ■'iv^jv+i - 2\\lN + 2' 

(2.25) 

The extension to d-dimensional systems was considered by Fiuràsek, Filip and 
Cerf [2005] who investigated the universal asymmetric quantum triplicator, which 
produces, from a single replica of a qudit, three clones with three different fidelities 
f^, and f'^'. A simple parametric description of the class of the optimal universal 
highly-asymmetric triplicators was provided, extending Eqs. (2.5) and (2.6). It was 
proved that the optimal fidelities can be expressed as 



1 - - 




N 



►1+1+1 



-1+1+1 
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d 

d-1 
d 



í3^ + 7' + 



2/^7 
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d-l 



■1+1+1 



2a(3 ' 



(2.26) 



where the positive real parameters a, j3, 7 satisfy the normalization condition 



+ + 7^ + + Q!7 + (3^) = 1. 



(2.27) 



2.10. CONTINUOUS- VARIABLE CLONING MACHINES 

Another interesting extension of quantum cloning, which is often termed 
continuous-variable quantum cloning, concerns the case of quantum systcms lying 
in an infinite-dimensional Hilbert space. In Cerf, Ipe and Rottenberg [2000], the 
cloning of the set of coherent states la) = e""*^ |0) was investigated, with |0) denot- 
ing the vacuuni statc, being the bosonic creation operator, and a = {x + ip)/V2 
being a c-numbcr which defines the position {x,p) oi \a) in phase space. Hcrc 
X and p are the so-called quadrature components. A set of l-to-2 (symmctric or 
asymmctric) cloning machines that are covariant with respect to the Wcyl group of 
displacements in phase space was derived. The symmetric l-to-2 Gaussian cloner 
was found to have the fidelity 

= 0.667 (2.28) 

and was conjectured to be optimal. It causes an independent Gaussian noise on x 
and p, with a variance equal to one shot-noise unit. Thus, the two clones are left in 
a thermal state (containing on average 1/2 thermal photon) which is displaced by 
a. Let us also mention the independent derivation of this l-to-2 Gaussian cloner 
as well as its extension to múltiple clones (M > 2) by Lindblad [2000]. 

Cerf and Iblisdir [2000] latcr on derived an upper bound on the fidelity of the 
symmetric A^-to-M Gaussian cloncrs, bascd on a link with state estimation theory. 
Since it coincided with Eq. (2.28) for A^ = 1 and M = 2, this proved that the 
above cloner is indeed the optimal cloner by means of a Gaussian operation. Cerf 
and Iblisdir [2001c] then showed that this l-to-2 Gaussian cloner can be realized 
simply by use of an optical parametric amplifier of gain 2, followed by a balanced 
beam splitter. The cloning noise then originates from the vacuum fluctuations of 
the ancillary modes that arc coupled to the input mode. 

In Cochranc, Ralph and Dolinska [2004], it was shown that if the cnscmblc of 
input coherent states has a finitc width, the l-to-2 Gaussian cloning can be achieved 
with a higher fidelity. Clcarly, if the task is to clonc a coherent state drawn from 
a distribution that is pcakcd around the origin of phase space, the vacuum state 
is a very good approximation of the original state, so cloning with a fidelity close 
to one is possible. The above fidelity fi^2 corresponds to the opposite situation 
of an infinitely wide input distribution, that is, an arbitrary input coherent state. 
If the input coherent state is distributed according to a Gaussian distribution of 
zero mcan and given variance, Cochrane, Ralph and Dolinska [2004] gave a closed 
formula for this fidelity as a function of the variance. 

The optimal A^-to-M Gaussian cloning transformation that achieves the above- 
mentioned upper bound was obtained in Braunstein, Cerf, Iblisdir, van Loock and 
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Massar [2001], Fiuràsek [2001a], yielding 



MN 



f^^-*^ - MN + M-N- ^^"^^^ 

As for discrete-dimensional states, these cloners tend, at the limit M oo, to the 
optimal measurement of with the fidelity 

/^loo = (2.30) 

The optical realization of these symmetric N-to-M cloners was also described there, 
while it was generaUzed to asymmetric l-to-2 cloners in Fiuràsek [2001a]. In the 
latter case, the balance between the fidelities of the two clones follows 

/a'' = T^' /b'' = T-^> cTAas = 1/2 (2.31) 

which corresponds to the no-cloning uncertainty relation derived in Cerf, Ipe and 
Rottcnbcrg [2000]. Hcrc, a\ and (T^ arc the variances of the added noise on clone 
A and B, while one shot-noise unit is taken as 1/2. 

Finally, in Cerf and Iblisdir [2001a], a more general class of Gaussian cloners was 
charactcrized, which transfornis rcplicas of an arbitrary coherent statc \a) and 
N' replicas of its phase-conjugate |q;*) into M clones of \a) and M' clones of |a*), 
with N — N' = M — M' . For well chosen ratios N' /N, this cloner was shown to 
pcrform bcttcr than the {N + N')-to-M cloner. In addition, the spccial casc of the 
balanced Gaussian cloner, with N = N' and M = M' was shown to be optimal 
among all cloners in this class in the sense that it yields the highest fidelity for 
fixed N + N' and M + M', namely 



I 

Interestingly, in the limit M — ^ oo, we have 



fS^N^M,M - 4^2^ ■ (2.32) 



/w3v^oo,oo — _|_ ^ (2.33) 

which means that the optimal measurement of |a)'^^|a*)'^''^ gives the same fidelity 
as the optimal measurement of \a)®'^^ , instead of |a)®^^ as a simple counting of 
states sccms to imply. This advantage of phase conjugation was first noted in Cerf 
and Iblisdir [2001b]. 

Coming back to the question of the symmetric N-to-M cloning of coherent states, 
Cerf, Krücgcr, Navcz, Wcrncr and Wolf [2005] havc rcccntly invcstigated the ques- 
tion of whether the above Gaussian cloners really provide the absolute highest 
fidelity or, instead, transformations outside the realm of Gaussian operations need 
to be considered. Against all intuitions, it was shown that, provided M is fi- 
nite, the cloning transformation that optimizes the single-clone fidelity is slightly 
non-Gaussian. For example, the optimal symmetric l-to-2 non-Gaussian cloner of 
coherent states was shown to have the fidelity 



(2.34) 
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strictly larger than f^^2 = 2/3 ~ 0.667. In contrast, the optimal cloners of coherent 
states with respect to the global fidelity remain Gaussian. This discrepancy between 
optimal cloners with respect to single-clone or global fidelities is reminiscent to the 
situation for phase-covariant cloners in finitc-dimcnsional spaces. 

For a review on continuous- variable quantum cloning, see Cerf [2003] and Braun- 
stein and van Loock [2005]. 



2.11. PROBABILISTIC CLONING MACHINES 

All the above listed cloning machines are deterministic, i.e., they always produce 
(imperfect) clones. However, one can also consider probabilistic cloning machines, 
which somctimcs fail to gcncratc the doncs but, if thcy succccd, thcn the cloncs 
exhibit higher fidelities than those achieved by the best deterministic cloners. The 
concept of probabilistic cloning was introduced by Duan and Guo [1998a], Duan 
and Guo [1998b], Chcflcs and Barnctt [1998], Chcflcs and Barnctt [1999] who in- 
vestigated cloning of a discrete finite set of pure states. They showed that a set 
of linearly independent states can be perfectly copied with some probability p. In 
particular, an cxact l-to-2 cloning of two generally non-orthogonal pure states 
and ['02) is possible with probability 

Pi^2 = ^^..] I , (2.35) 

The probabilistic cloning was then extended to infinite continuous sets of input 
states in Fiuràsek [2004]. It was shown that the optimal universal cloning cannot 
be improved by using a probabilistic cloning strategy, due to a vcry high underlying 
symmetry of the problem. Nevertheless, if one considers cloning of some restricted 
set of states, then probabilistic cloning may become useful. A particular exam- 
ple is the optimal N-to-M phase-covariant cloning of qubits, where the optimal 
probabilistic cloner achieves the single-clone fidelity 



/.pc, prob/QX _ 1 



M 



k=0 



k+[\{M-N)\ 



whcrc [x\ denotes the intcgcr part of X. ¥ovN>l the fidelity /^!:,5^°^(2) is larger 
than the fidelity ff^_^]^{2) of the optimal deterministic phase-covariant cloning. 



2.12. ECONOMICAL CLONING MACHINES 

The l-to-2 cloning transformation for rf-dimensional systems (qudits) can typically 

be expressed as a unitary opcration on the Hilbert spacc of thrce qudits — the 
input, a blank copy, and an ancilla. The presence of an ancilla significantly aífects 
the experimental implementation of the cloning operation, which becomes more 
complicated and sensitive to decoherence. Thesc problems, which might drastically 
reduce the achieved cloning fidelity, may significantly be suppressed if an "econom- 
ical" approach is followed, which avoids the ancilla. The l-to-2 cloning is then 
realized as a unitary operation on two qudits only: the input and the blank copy. 
This is obviously much simpler to implement because it requires less qudits and 
two-qudit gates, and it requires to control the entanglement of a pair of qudits 
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only. It is thus likcly to be much less sensitive to noise and decoherence than its 
three-qudit counterpart. 

To date, the only l-to-2 cloning machine for which an economical realization is 
known is the phasc-covariant qubit cloner duc to Niu and GrifRths [1999], which 
optimally clones all statcs on the equator of the Bloch sphere, = 2~-^/^(|0)^ + 
e"^|l)^). The qubit to be cloned is coupled to another qubit which becomes the 
second copy and is initially prepared in state |0)b. The unitary two-qubit trans- 
formation reads 



for a symmetric phase-covariant cloner. It can easily be extended to an asymmetric 
setting as wc will show in Eq. (6.4). 

The possibility of economical realization of various l-to-2 cloning machines for 
qudits has been analyzed in detall by Durt, Fiuràsek and Cerf [2005]. It was shown 
therc that economical universal cloning is not possible for any d. It was also argued 
that the optimal l-to-2 phase-covariant cloning of qudits docs not admit economical 
implementation for any d > 2, and this assertion was rigorously proved for d <7. A 
suboptimal economical phasc-covariant cloner was nevertheless constructed, which 
does not require an ancilla and achieves the fidelity 



which is only slightly below that of the optimal cloner. Similarly, it was argued 
that the l-to-2 Fourier-covariant cloning cannot be realized economically, albeit in 
dimension d = 2 (in which case it is unitarily equivalent to the phase-covariant 
cloner). 

The concept of economical cloning can be extended to A^-to-M machines. As 

dctailcd in Scction 6.1, the paper by Fan, Matsumoto, Wang and Wadati [2001] 
implies that the optimal N-to-M phase-covariant cloning of qubits {d = 2), which 
maximizes the single-clone fidelity, admits an economical implementation for any 
N and M > N. Morcovcr. the economical phasc-covariant cloning of d-dimcnsional 
systems (qudits) is also possible provided that M = kd + N, where k is an integer, 
see Buscemi, D'Ariano and Macchiavello [2005]. 

§ 3. One-to-two quantum cloning as a CP map 

3.1. ISOMORPHISM BETWEEN CP MAPS AND OPERATORS 

A very useful characterization of cloning relies on the isomorphism between com- 
pletely positive maps S : Tíin — > Tíout and positive semidcfinite operators S' > 
acting on Hhi $5 Wout, whcre Hin and Tíout denote, respectively, the input and out- 
put Hilbert spaces of S, see Jamiolkowski [1972], Choi [1975]. To construct this 
isomorphism, consider a maximally entangled state on H^, 



\0)a\0)b ^ |0)a|0)s 



|i)a|o)s-4(|o)^|i)b + |i)a|o)b), 



(2.36) 




(3.1) 
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where d — dim(7íin). If the map S is applied to the second siibsystcm of 1$+) while 
the first one is left unchanged, then the resulting (generally mixed) quantum state 
is isomorphic to <S and reads 

S=ÇI(^S)d^+, (3.2) 

where = |$+)(<J>+| and X stands for the identity map, while the prefactor d is 
introduccd for normahzation purposes. The map <S can be characterized in terms 
of the state S as follows, 

Pout = 5(pin) = Tri^iipïn Iont)S], (3.3) 

where "in" labcls the input spacc, / is tlic identity opcrator, while T denotes the 
transposition in the computational basis. A trace preserving map S implies that S 
satisfies the condition 

Trout[5] = /in, (3.4) 

while the completo positivity condition on S translatcs into S > 0. 

In the foUowing, we shall make this description specific to the one-to-two quantum 
cloning machines, which produce two copies of a single rf-dimensional system (often 
called a qudit), sec Fiuràsck [2001b]. The output Hilbcrt space is endowcd with a 
tensor product structure, Hout = 'Ha '^'Hb, where the subscripts A and B label 
the two clones. For each particular input state \tjj), the joint state of the clones is 

5(V) =Trin[(V'í;®/^B)^], (3.5) 

where ip = \ip){tp\. It will be useful in the foUowing to note that V'í^ is a rank- 
one projector onto the state IV'in)' where * denotes the complex conjugation in the 

computational basis. 

Using Eq. (3.5), the fidclity of the clones A and B is givcn by 

Fa{S,^) = Tt[{,Pa®Ib)S{iP)] = Tt[Í^I®^a®Ib)S], 

FB{S,tP) = Tr[(/^®^B)5(^)] = Tr[{tl;l(g>lA'E>^PB)S]. (3.6) 

The symmetric cloning machines are defined as the maps iS verifying FA{S,tp) = 
Fb{S,4>), VV». Otherwise, the cloning machines are called asymmetric. When con- 
sidcring a universal cloning machine, we require that both Fa{S,i})) and FbÍS^iP) 
are independent of ■0, for all state ■0 in Tíin- Other cloning machines, such as the 
phase-covariant, Fourier-covariant, or real cloning machines will correspond to a 
constant fidelity over a restricted set R of input states 

3.2. COVARIANCE CONDITION 

In what follows, we will always assume that the set R of input states ij} is invariant 
under the action of the group G(fi) of unitàries {Í7a,|a; e O}, that is 

U^RUl = R, Vwefi. (3.7) 

The universal cloning machine is the special case Ua, G SU(cí). A useful figure of 
mèrit to measure the quality of cloning is the global fidelity. 



F{S,i>) = Tr[(VA«)VB)5(V)] = Tr[{i,l^i,A^Í'B)S], (3.8) 
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which measures how wcU thc joint state of the two clones approximates í/;**^. When 
looking for a cloning machine that optimally clones all the states of set R, one 
generally defines the cloning fidelity of map <S as the infimum of the global fidelity 
over all input states ip, 

F{S) = inf (3.9) 

It has bccn shown by Werncr [1998] that. by using thc so-callcd twirling operation, 
there is no loss of generality in assuming the optimal cloning machine to be covariant 
with respect to the group G{Q), hence the cloning fidelity to be state- independent 
within thc set R. Thc twirling operation consists in randomly applying a unitary 
Uu to the input state and then undoing this by applying the reverse unitary to 
each of the two clones with the probability density dco equal to the Haar measure 
on the group G{Cl). This results in the twirled map 

M = / S^[i;]du (3.10) 
Jn 

with the rotated map S^^ being defined as 

Su^m = Ul^^S [Uu^^Ul] UT- (3.11) 
The core of the argument is that 

F{S) = inf F(5,^) < inf / F{S,U^4>Ul) dw 

= iuí í F {S^, (/.) dw = inf í^(5twiri, (p) = í^(5twiri), (3.12) 

where we have uscd thc invariancc of R under the unitàries t/^^, the invariance of 
the trace function under and the linearity of the fidelity in S. As a result, the 
operation of twirling can only increase the cloning fidelity, so that the twirled map 
5twiri is at Icast as good as cach of its constituent maps S^ü ■ Finally, as mcntioned 
earlier, we note that »Stwiri is covariant with respect to the group G(íi), that is 

St^iAMul] = u^^St^iri{ml^^ € e R. (3.13) 

Physically, this covariance property means that rotating the original state is exactly 
equivalent to rotating the two clones by the same amount. This also implies that 
F{Stwiri,^) does not depend on ip. 

In summary, we have shown that when looking for an optimal quantum cloning 
machine (i.e., a machine that maximizes the worst-case global fidelity), it is suffi- 
cient to consider cloning maps that are covariant with respect to the group under 
which the set of input states R is invariant. The cloning fidelity is therefore state- 
indcpcndcnt within the set R. In Keyl and Werner [1999], it was proven that this 
reasoning also applies more generally to the quantum cloning machines that maxi- 
mize the single-clone fidelities {Fa and Fg) instead of the global fidelity, provided 
that universal cloning machines are considered. This, howcvcr, docs not hold for 
all quantum cloners (see, e.g., the case of phase-covariant or continuous- variable 
cloners treated later on). 
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3.3. CLONING AS A SEMIDEFINITE PROGRAMMING PROBLEM 

Corning back to the characterization of the map S via its associatcd operator S, 
we can now use the fact that the optimal cloning machine must have a state- 
independent fideUty over the considered set of input states. One can then turn to 
the average performance of the cloning machines, which is measured by the mean 
fidelities, 

F^(5)= / Fb{S)= i FB{S,^)diP, (3.14) 

where the measure determines the kind of cloning machines we are dealing with. 
In particular, universal cloning machines correspond to choosing dijj to be the in- 
variant measure on the factor spacc SU (d) / SU{d— 1) induced by the Haar measure 
on the group SU{d). The fidelities (3.14) can be expressed as linear functions of 
the operator S, 

Fa = Tr[SRA], Fb = Tr[SRB], (3.15) 
where we have defined the positive semidefinite operators 

Ra= il^l,®il)A®lBdil}, Rb= Íjfn<»lA<»Í}Bd^. (3.16) 

We will see in the next Sections how these operators can be calculated for different 
kinds of cloning machines. Note that for a symmetric cloning machine, one should 
simply maximize the average of the mean fidelities, 

= 1[Fa{S) + Fb{S)] = Tv[SR], (3.17) 

with R = {Ra + Rb)/2. This can be justified by an argument similar to the one 

used for the twirling operation. By averaging over the permutation b(rt.w(H^n the 
two clones, one obtains a map whose global fidelity can only be better than that of 
the original map. Therefore, we can restrict ourselves to cloning transformations 
that are covariant with respect to the interchange of the clones, hence satisfying 
Fa{S) = Fb{S). 

From this argument based on twirling and permutation, wc conclude that maxi- 
mizing the mean fidelity, averaged over the two clones (with cqual wcights) , should 
yield a cloning map which has a state-independent and clone-independent fidelity. 
The asymmetric cloners can also be obtaincd with the same maximization but by 
putting different weights in front of Fa and Fb- 

An interesting point to note is that finding the optimal cloning map S reduces 
to a semidefinite programming problem, namely finding the operator S verifying 
6* > and TrAsíS*] = /in that maximizes Tr[SR], with R depending on the consid- 
ered cloning machine (Audenaert and De Moor [2002]). Therc cxist very efíicient 
numcrical mcthods for solving semidefinite programs, sec c.g. Vandcnbcrghe and 
Boyd [1996]. Evcn morc importantly, it can be shown with the liclp of Lagrange 
duality lemma that the optimal cloning trace-preserving CP map, which maximizes 
Tr[5iï], must satisfy 

(fí - Ai„ ® /ab)5 = 0, (3.18) 
Ain ® /ab - i? > 0, (3.19) 
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where A > is a positive semidefinite operator whose matrix elements represent 
the Lagrange multipliers accounting for the trace-preservation constraint Tr 43 [S] = 
Jin. Note that A can be expressed in terms of the optimal CP map, A = Ttj\b[S R]. 
If both equations (3.18) and (3.19) arc satisfied, thcn S is the optimal CP map 
maximizing Tr[S'iï], a property which is uscful to prové and check the optimality 
of a givcn map S that is conjectured to be optimal. 

The proof that Eqs. (3.18) and (3.19) imply optimality is rather simple and we 
briefly sketch it here. Suppose that (3.19) is satisfied, then it holds for any trace- 
preserving CP map that Tr[S'(A ® / - i?,)] > and Tr[A (g) IS] = TrA, due to the 
trace-preservation condition. It follows that the fidclity is \ipper bounded by the 
trace of the Lagrange multiplier, Tr[RS] < Tr[A], and that the optimal map which 
satisfies (3.18) satinates this bound. 

Very often a simpler method is sufRcient to prové the optimality, namely, the 
fidelity can be bounded by the maximum eigenvalue r^ax of R. Since R < r^^J, 
we immediately have 

F{S) < d (3.20) 

where d = dim(7íin)· If there exists a CP map S which saturates (3.20), then this 
transformation is optimal. Note, however, that in certain cases such as the cloning 
of a pair of orthogonal qubits, the bound on fidelity (3.20) is not tight and cannot 
be saturated. 

3.4. DOUBLE-BELL ANSATZ 

Let us now consider the unitary realization of the cloning map S. We know that 
any CP map can be realized physically by supplementing the input systcm with an 
ancilla (hence, extending the Hilbert space) and acting with a unitary operator in 
this extended space. Here, the ancilla can be viewed as the cloning machine itself, 
and it must be traced over after applying the unitary operator. The resulting map 
can be written as 

5(^i„) = Tri„,c[(V',^ ® lABc)n, (3-21) 

where C denotes the cloning machine and E is the operator that is isomorphic 
to this extended map Tijn TLa ® "Hb ® 'Hc ■ Since this extended map is some 
unitary operation Uabc in the extended space Ha ^Hb® Hc, the operator S 
must be some (imnormalizcd) rank-one projector or pure state in the joint space 
Hin ® Ha ® (g) Tic- We thus have 

S = |cr)(cr|, with |a·)in,ASC = Vd (7in ® Uabc) |$+)in,A|0)s,c (3.22) 

where \Çi)b,c is the (arbitrary) initial state of the blank copy B and cloning machine 
C, and |<ï>+) is defincd as in Eq. (3.1). The prefactor with d = dim(7Yin), is in- 
troduced for normalization purposes, and the cloning map corresponds to Eq. (3.5) 
with S = TrcS. 

Physically, the state \a) has a very simple interpretation, sec Fig. 2. If we start 
with two qudits prepared in a maximally entangled state |$^) and process one of 
them in the quantum cloning machine while the other one is left unchanged (kept 
as a reference), then \a) is the joint state of this "reference" qudit (denoted as "in" 
since it keeps a memory of the input state), the clones A and B, as well as the 
cloning machine C. Remember that if we project the reference qudit onto the state 
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Figure 2: Preparation of the state a fully characterizing the cloning transformation. 

The input of the cloning niachine is maximally cntanglcd with a rcfcrcncc qudit 
labeled "in" . The two clones are contained in the outputs A and B while C refers 
to an ancilla of the cloning machine. 



li/"*), then. in the abscncc of cloning, qudit A is found in state \ip). By causality, 
it is irrelevant whether this projection onto is done before or after the cloning 
machine has been applied on qudit A. Therefore, projecting the reference qudit 
of state |(t) onto \·4)*) yields the joint state of A, B, and C that would have been 
obtained by cloning the state jf/'), namely 

IV') \Í>out) ABC = in(V'*|o·)in,ASC· (3.23) 

We can say that \a) fully encodes the Information about the cloning of any state. 

It was suggested in Cerf [1998], Cerf [2000a], Cerf [2000b] that a genèric form for 
state involving a superposition of double-Bell states may encompass most of the 
interesting quantum cloning machines, including the universal or state-dcpendent 
- symmetric as well as asymmetric - cloners. This so-called double-Bell ansatz 
corresponds to taking \cr) = \/d\A), with 

d-l 

\·A)in,A;B,C = "^n,n \^m,n) in,A\^m,n) B,C , (3.24) 

m.n— 

where it is asstuned that it is sufíicient to use a Hilbert space for the cloning machine 
C which has the same dimension d as that of the input or the clones. Note that 
the Schmidt decomposition of |^) for the partition "in" vs. ABC implies that 
dïm{Tíc) > d- In Eq. (3.24), the a^ n's are complex amplitudes which satisfy the 
normalization condition X^m n=o l'^ni.nP = 1) while \^m.n) denote Bell states in d 
dimensions. As the latter states play an important role in what foUows, we will 
first discuss them in detalls, as well as some useful related properties. Note also 
that 1^) is a quantum state of norm one, while \a) has a norm equal to d. 

3.4.1. Useful properties of rf-dimensional Bell states 

A Standard generalization of the Bell states in d dimensions is 

\^m,n) = ^^7"^|j)ib + m>2, (3.25) 
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whcrc 1 and 2 denote two d-dimensional systems. Thcse states form a set of (P 
maxiïnally-entanglcd states of systems 1 and 2, where m,n Çi {0, l, ■ ■ ■ d — 1} and 
7 = £;27ri/d s^ands for the d-th root of unity. Note that, in what foUows, the "bra" 
and "ket" labcls arc always taken modulo d. In the case of qubits {d = 2), we 
recover the Standard Bell states 

|$o,o> = |$+) = (|00) + |11»/V2, Ic&oa) = 1$-) = (|00> - |ll»/%/2, 
|$i,o) = 1*+) = (|01) + |10))/V2, = I*-) = (|01) - |10))/V2. 

(3.26) 

Taking the partial trace of any state í'm.n = |^m,ra)(^m,rí| over one of the two 
systems (1 or 2) results in the maximally mixed state, 

Tri($™,„) = Tr2($™,„) = I/d, Mm, n (3.27) 

so that the states |$m,n) are indeed maximally entangled. It is easy to check that the 
states \^m,n) form a complete orthonormal basis in the considered rf^-dimensional 
Hilbert space. The resolution of identity reads 

m,n=0 m,n jj' 

= Y.Y. \j){f\ ® li + m){j' +m\= h2, (3.28) 

"» 3,j' 

where we have used the identity ^ X^^=o T""' ~ ^ófi- 
Let us focus on the Bell state with m = n = 0, that is 

l*o.o) = ^Eb')b') (3.29) 

which is another notation for the state |<3>"'") as defined in Eq. (3.1). This state is 
particularly useful because it satisfies the relation 

(i7*0[/)|$o,o) = l^o.o) (3.30) 

for any unitary transformation U, as can readily be checked by using the unitarity 
condition UW = I and the completeness relation \ = I- Note that the 
symbol * denotes the complex conjugation operation in the computational basis 
{|j)}; thus |j*) = The identity (3.30), or equivalently 

(/® [/)|$o,o) = (C/^®/)|$o,o) (3.31) 

corresponds to the foUowing useful property: if the joint system 12 is prepared in 
the state [í'o.o) and system 1 is projected onto \^p*), then the resulting state of 
system 2 is IV»)- Indeed, taking IV») = í^|0), we have ('^*| = (0|{7'^, so that 

{\r)m ® /)i*o,o) = (i^*)(oi ® i){u^ ® /)i$o,o) 

= (|V*)(0|®/)(/®[/)|$o,o)=rf-'/'|V'*)(t^|0)) = íí-^/'|V*)l^)· (3.32) 
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Interestingly, it makes no diffcrcncc whcther 1 prepares and scnds the state jí/i) to 
2, or whether 1 projects its part of a shared entangled state \^o,o) onto \tp*) so to 
create \ip) at a distance on 2. 

In what foUows, we will also nccd thc discreta group of Weyl-Heisenberg operators 
(also called error operators), namely 

d-l 

í^™,n=^7^"IÍ + W01 (3-33) 
j=o 

with m,n G {0,1, •••d— 1}, generalizing the Pauli matrices for morè than two 
dimensions. For qubits {d = 2), we have 

Eo,o = I, -2^0,1 = o-z, 

Elfi = Cx, -E'l,! = (Txf^z = —icTy. 

(3.34) 

In arbitrary dimension, the error operator Em^n shifts the state by m units (modulo 
d) in the computational basis and multiphes it by a phase so as to shift its Fourier 
transform by n units (modulo d). Indeed, in the computational basis {\j)}, we have 
Em,o\j) = \j + m), while in the dual basis {\\j)) — {í/Vd)J2k=o^^''\^)}^ ha.ye 
Eo,n\\j)) = lli + n)). The error operators fulfiU the foUowing properties, 

Em,n = Em-n, E^^^ = 7"™" (3.35) 

Em,n = 7™"-2'-m,-n5 Em^nE^^u = 7"'*-£'rn+/n,n+i/· (3.36) 

Interestingly, the Bell states can be transformed into each other by applying an error 
operator locally (on one of the two systems, leaving the other one unchanged), 

\^m,n) = (/ Em,n)\^0,0) = {E^,n 8» ^)I^O,o). (3.37) 

This also implies that the Bell states are invariant (up to a phase) under correlated 
error operators {E* ^, (g) E^^^). We can check this by calculating 

{El^^ (g) E^^^) |$„,„) = {I® £;^,^S™,„) {El^^ (g) /) |$o,o) 

= (j (g) í;^,.í;„,„) (/ ® El,) |$o,o) 

= 7— ""(/ £„,„) |$o,o) = 7™"""'' l*m,n), (3.38) 
where we have used property (3.37) as well as 

7? c; z?t _ ^fíf E E E — ^t^'^-'^t^ E E 

3.5. HEISENBERG CLONING MACHINES 

Coming back to the double-Bell ansatz, Eq. (3.24), the quantum cloning machine 
is thus completely characterized by the dx d matrix a = {am,n}· The form (3.24) 
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is particularly interesting because, when tracing over B and C, the systems "in" 
and A are left in a mixed state that is diagonal in the Bell basis, 

d-1 

Pin,A = ^ \am,nf ^m,n (3.40) 
m,n=0 

with $m,ri = \^m,n){^m,n\- Since the Original system is maximally entangled with 
the reference system "in" (the initial state being |$o,o)), this implies that clone A 
undergoes the error Em,n with probability Ioto,?»!^- It emerges from a Heisenberg 
channel characterized by the probability distribution |a„j,„p. 

An important property of state \A) is that, when interchanging clones A and B, 
it can be re-expressed as a superposition of double-Bell states albeit with different 
amplitudes, 

d-1 

\A)in,B;A,C = Yl W I *m,n)m,B | A,C (3.41) 

m,n=0 

with 

^ d-1 

bm,n = 2 12 ^"""'"'«x,,- (3.42) 
x,y=0 

Again, when tracing over A and C, systems "in" and B are left in a Bell-diagonal 
mixed state, 

d-1 
m,n=0 

implying that clone B undergoes the error Em,n with probability |òm,„P (it emerges 
from another Heisenberg channel). Remarkably, Equation (3.42) implies that the 
matrix b = {6m,ri} is rclated to a = {am,n} by a (bivariate and d-dimensional) 
discrete Fourier transform, b = J^[a\. So, the cloning map can be characterized 
equivalently by the matrix a (characterizing the noise of clone A) or its Fourier 
transform b (characterizing the noise of clone B), and we see that the comple- 
mentarity between these two clones simply originates from a Fourier transform; 
the more noisy is clone A, the less noisy is clone B. This leads to a no-cloning 
imccrtainty relation, sec Ccrf [1999], Ccrf [2GGGb]. 

Finally, we can use the ansatz (3.24) to express the map associated with an 
arbitrary Heisenberg cloner in the simple form 

d-1 

IV") IV'out) = Y "'"."Trin 

m,n=0 

d-1 

= am,nEm,n\'4^)A<»\^*m,n)BC- (3.44) 

m,n=0 

Incidentally, we note here that by measuring the clone B together with the cloning 
machine C in the Bell basis, we get a pair of indices (m, n) which can be used to 



(ipfn ® IaBc) ({^m,n)in,A ® (C„)bc) 
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undo the noise on clone A simply by applying E^^ This process, which bears 
some analogy with quantum teleportation, will be exploited in Section 5.3 in order 
to convert a symmetric cloner into an asymmetric cloner. 

3.5.1. Covariancc w.r.t. thc Wcyl-Hciscnbcrg group 

It can be proven that the Heisenberg cloning machines are covariant with respect 
to the discrete Weyl-Heisenberg group of error operators {Em^n}- Recali that Ej^^ 
corrcsponds to a cyclic rclabcling of thc computational basis states, while Eq ,^ 
corresponds to a cychc relabehng of the dual basis states; Ej^^^ = E^^qEq^^ sim- 
ply corresponds to a sequence of these cyclic permutations. Thus, the Heisenberg 
cloncrs arc covariant with respect to cyclic permutations of the basis states in the 
computational and dual basis. 

Using Eq. (3.5), it can be easily shown that the covariance condition of the 
cloning map <S with respect to the unitary operator U, namely 

5[Í/V'í7t] = [/®25(^);7t®2 VV» e (3.45) 
translates into the condition 

{U* Í7®2) g ([|T ^ jjt®2) ^s^[S,U*® U®^] = (3.46) 

on thc operator S that is isomorphic to S. Wc may also impose that when the 
original is transformed according to the unitary U, the cloning machine is trans- 
formcd according to thc unitary U* . This condition, named strong covariance, can 
be expressed as a constraint on state \a) or |^), namely 

(U* (g> U^^ (g> U*)\A)in,A,B,C = \A)in,A,B,C- (3.47) 

It was shown recently that, provided that the set of input states is invariant with 
respect to the Weyl-Heisenberg group, the class of strongly-covariant cloning maps 
is equivalent to thc class of extremal covariant maps, sec Chiribclla, D'Ariano, 
Pcrinotti and Cerf [2005]. Thus, substituting covariancc with strong covariance 
grcatly simplifies the search for optimal cloners since, given that the covariant 
cloncrs form a convex set, it is sufEcient to search among extremal cloners. 

Thc strong covariance of the Heisenberg cloners can be checked by using condition 
(3.47) with U = E^^v for all ^ and v. This equation indeed holds for each component 
of A, namely 

= \^m,nWm,n), (3-48) 

where we have used Eq. (3.38). Thus, the Heisenberg cloning machines defined 

by the ansatz state |^) for an arbitrary matrix a havc thc nice property that they 
keep the same form when making a cyclic permutation of thc basis states (in both 
the computational and dual bases). 

This covariance property also implies that the reduced cloning maps are unital. 
It is trivial to prové that applying an error operator Em,n chosen at random (uni- 
formly among the (P possibilities) on an arbitrary state p always gives a maximally 
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disordered state, 

1 '^'^ I 

^ J2 Em,npEl^^^-. (3.49) 

m,n=0 

Consider an arbitrary input state of the cloner tp, the covariance and the linearity 
of the reduced cloning map <S^ or Sb imply that 

/ ^ d-i \ ^ d-i 

'^^'^ b2 E Em,nÍ^El^n\=-^ E E^,^ Sa,bW El,r^ (3-50) 
\ m,n=0 / m,n=0 

SO that, using Eq. (3.49), we verify that the reduced cloning maps «S^,b are indeed 
unital 

SA,B[I/d\ = I/d. (3.51) 



3.6. THREE SPECIAL CASES OF HEISENBERG CLONERS 
3.6.1. Universal cloners 

Let us now discuss several interesting special cases of Heisenberg cloning machines. 

The first example is the universal (or isotropic) cloning machine, where the channel 
underlying each output is a quantum depolarizing channel. This implies that all 



the probabilities p„ 



except po,o must be equal. The same holds for the 



probabilities qm.n = l&m.nP associated with the second clone. These conditions put 
very strong constraints on the matrix a, whose elements can be thus parametrized 
by two real coefficients v and x, 



(u - x) 6n,0 (5m,0 + X 



( V X 

X X 



X \ 

X 



\ X X ■ ■ ■ X j 

The Fourier transform yields the matrix elements of ò, namely 
with 



V — X 



X = 



{v' - x')Sn,oSm,0 + x', 

v + {cP - l)x 



v' = 



d ' d 
The cloning is a trace-preserving operation so the condition Tr[pin,A] 
satisfied, which provides the normalization constraint. 



l)x' = 1. 



(3.52) 



(3.53) 

(3.54) 
1 must be 

(3.55) 



If the input state is |0), the error operators Eo^n, Vn, leave it unchanged up to a 
phase, while all the other Em,nS produce a state that is orthogonal to it. Therefore, 
the fidelities of the two clones for any input state can be expressed as 



FA=v^ + {d- l)a;^ Fb = v'^ + {d- l)x' 



„'2 



(3.56) 
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Note tliat there is only a singlc free parameter x which controls the asymmetry of 
the cloner. We can also characterize the cloner as a function of the íidehty Fa oi 
the first clone, namely, 



1-Fa 



{d+l)FA-l 



d{d-iy d 
The symmetric cloner is obtained by putting x = x' , which results in 



1 



2d{d+iy 



d+1 
2d 



(3.57) 



(3.58) 



and is associated with the fidclity givcn in Eq. (2.3) 

Note that, as rigorously proved recently for any d in Fiuràsek, Filip and Cerf 
[2005], this isotropic Heisenberg cloner represents the optimal asymmetric cloning 
machinc which, for a fixcd fidclity Fa of the first clonc, maximizes the fidclity of 
the second clone Fb- Note also that the optimality of the Heisenberg cloners, based 
the double-Bell ansatz, was explained in Chiribella, D'Ariano, Perinotti and Cerf 
[2005] as a conscqucncc of the cxtrcmality of thcsc cloners. It is worth strcssing 
that by exploiting this double-Bell ansatz, these machines can be derived almost 
without any effort as they foUow immediately from the general isotropy and trace 
preservation conditions. 



3.6.2. Fourier- covar iant cloners 

As a second example, we shall consider the Fourier- covariant machine, which clones 
cqnally well two mutually unbiased bases, the computational basis {\k)} and the 
dual basis 

= ^ye2^'('='/'')|fc). (3.59) 
vd 



fe=o 



The cloner copies equally well the states of both bases if the matrix a has the form. 



= {v -2x + y) Sm,o Sn,o + (x - y) {6m.,o + Sn,o) + y 
( V X ■ ■■ x \ 

X y ■■■ y 



\ X y ■■■ y J 
where the parameters satisfy the trace preservation condition 

+ 2(d - l)a:2 + (d - 1) V = 1 
The matrix 6 then has a similar form with v, x, and y being replaced by 

v' = ^[v + 2{d-l)x + {d-iyy], 

x' = ^[v + {d-2)x + {l-d)y], 



(3.60) 



(3.61) 



(3.62) 



y' = 2^v-2x + y]. 
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The fidclitics of thc two clones are again given by Eq. (3.56) but, herc, wc have two 
free paramctcrs, say, x and y. To eliminate one of them, one has to maximize Bob's 
fideUty Fb for a given value of Ahce's fidehty Fa (using the normahzation rela- 
tion), which is a simple constrained optimization problem. The resulting optimal 
asymmetric cloner is characterized by 



Fa{1-Fa) i -Fa 

v = Fa, x=^ , y=^^. (3.63) 

which depends on the single parameter Fa ■ The symmetric Fourier-covariant cloner 
can again be obtained by setting x = x' and y = y' , which gives Eq. (2.16) for the 
fidelity. 



3.6.3. Phasc-covariant cloners 

As a third example, consider the phase-covariant machine, which optimally clones 
all balanced superpositions of the form 

= d-^/^[\0) + e^-^i |1) + • • • e"''-'-' \d - 1)] (3.64) 

where the (j)i's are arbitrary phases. Here, it can be easily shown that the matrix a 
must take the form 



am,n = (v-y) S,n,0 Snfi + {y - x) Sm,0 + X. 

/ V y ■■■ y \ 



\ X ■ ■ ■ X j 

while the trace preservation condition is 

+ (d - l)y2 + did - \)x^ = 1 
The matrix h has the same form, albeit with a;, and y being replaced by 

v' = ^[v + d{d - l)x + {d - l)y], 



(3.65) 



(3.66) 



(3.67) 



y 



-[v - dx + {d - l)y]. 
d 



If the input state is 1 10)) (all the states of the dual basis are balanced superpositions), 
the error operators Em.o, Vm, leavc it unchangcd up to a phase, while all the other 
Em.nS produce a state that is orthogonal to it. Thcrcforc, the fidelities of thc two 
clones are again given by Eq. (3.56), and we have two free parametcrs, say, x and 
y. We can eliminate one of thcm by maximizing Fb for a given Fa, which yields 
the optimal phase-covariant cloner. In the special case of a symmetric cloner, we 
have X = x' and y = y' resulting in the fidelity given by Eq. (2.17). 
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§ 4. N-to-M universal qucintum cloning 
4.1. OPTIMAL CLONING TRANSFORMATION 

In this Section we will focus on the universal (state-independent) cloning. An ideal 
universal N ^ M quantum cloning machine would be a device that prepares M 
exact clones of an arbitrary state tp £ H from N copies of tp- The input Hilbert 
space of the cloning transformation is the symmetric subspacc of N qudits, 

and d = dimH shall denote the dimension of the Hilbert space of the input states 
in what follows. As already explained above, exact deterministic quantum cloning 
is forbidden by the linearity of quantum mechanics and only approximate copying 
with fidelity less than unity is possible. 

As noted before, two different kinds of cloning fidelities are considered in the 
literature. The global fidelity compares the global state of M clones with the ideal 
output tp'^^ . Let S denote the cloning CP map. Then the global fidelity of cloning 
of the state V can be expressed as F^^'^^iS, ^) = Tr[V'®"<S(V'®^)]. Generally, the 
fidelity of the cloning can dcpend on íp and one may define the cloning fidelity as 
the infimum of F^_^^{S, %p) over all input states ip, 

FT^^{S) = inf TV[V®^5(V^)]. (4.1) 

The single-clone fidelity quantifies how well each clone resembles the desired output 
iP. For A;-th clone we can write F^''_!^]^'^{S,^,k) = TrliPTr^S {^'^^ )] where Tr'i, 
denotes trace over all M qudits except for the k-th qudit. When judging the 
performance of the cloning machine, we should take the infimum of F^_^j^'^{S, ip, k) 
over all input states and all M clones and define 

FT^i'^iS) = inf inf Tr[VTrí,5(V^)], (4.2) 

k yj 

where k e {l,---,!^}. 

The universal cloning machine should clone equally well all quantum states so 
the fidelity should not depend on ip. Any transformation S can be converted into 
a universal cloning transformation whose fidelity is state independent by twirling 
operation consists of applying randomly a unitary U (O) to each input ip and then 
undoing this by applying a unitary W{íl) to each clone with the probability density 
díl equal to the Haar measure on SU(2), see also Sec. 3.2. The effective map 

5twiri(V')= / C/^®'^(í2)5[(C/(íl)VÍ7t(n))®^][/®^(n)dí2 (4.3) 
Jn 

is covariant, i.c, St^i,i[(UtpU^)'^'^] = U'^^ St^iri{i^'^^)W^'^' and, consequently, 
Fn^^m does not depend on -0. To guarantee the independence of the single-clone 
fidelity on the clone index k, it is also necessary to randomly permute the M 
clones after the twirling. The important feature of the twirling operation and the 
pcrmutations is that they do not modify the mean fidelity calculated as the average 
of F^^^í^iS, V) or ^ E^i Fril^^'iS, ^, k) over all input states |V) = t/(í2)|Vo) 
with the measure díi. 

The universal cloning has been studied extensively by many authors (Buzek and 
Hillery [1996], Gisin and Massar [1997], Hillery and Buzek [1997], Buzek and Hillery 
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[1998a], Buzek, Hillery and Knight [1998], Ccrf [1998], Werner [1998], Niu and 
Griffiths [1998], Keyl and Werner [1999], Ccrf [1999], Ccrf [2000a], Ccrf [2000b]). 
The task of cloning can be rcphrased as diluting thc quantum information carried 
by the N input qudits into M output qudits. The universal cloning should not 
prefer any direction in the Hilbert space and should be isotropic. As shown by 
Werner [1998], the optimal universal cloning operation <Sopt can be expressed as 
foUows, 

= ^(S^"^.'^^^^'' ® ^^^^-^^)nl.,„ (4.4) 

where ^ is the projector onto the fully symmetric (Bose) subspace of M qudits 
and 

D(M,,)^C^^-') (4.5) 

is the dimcnsion of this subspace. We can see from Eq. (4.4) that the optimal 
cloning formally consists in attaching M — N blank copies prepared in maximally 
mixed state I/d to the input state r/»®^ and then projecting the whole state of M 
qudits onto the symmetric subspace of M qudits, Tíf^. With proper normalization 
as given in Eq. (4.4) iSopt is a trace preserving completely positive map and can be 
therefore realized deterministically. 

The maximal global cloning fidelity achieved by the optimal cloner (4.4) reads 



puniv,G ^ DiN,d) ^ {d + N-l)\M\ 

D{M,d) {d + M-iy.N\' ^ ^ 



The density matrix of each output clone is a convex mixture of the input state tp 
and the maximally mixed state I/d, 

p = j^+l{l-j)L (4.7) 

This cxprcssion rcvcals thc high isotropy of thc universal symmetric quantum 
cloning which is fully characterized by a single parameter, namely the shrinking 
factor 7(7V,M), 

,(iV,M) = -^^. (4.8) 

^ ^ N+d M ^ ' 

The single-clone fidelity can be immediately dctcrmincd from Eq. (4.7) and we get 

,univ,SC _ MjV + M + jV(d -l) 

M(iV + d) 



4.1.1. Conncction with quantum state estimation 

There is a close relationship between optimal quantum cloning and optimal quan- 
tum state estimation. As shown by Bruss, Ekert and Macchiavello [1998], in the 
limit of infinite number of cloncs, M oo, the singlc-clonc fidelity F^_^j^'^ bc- 
comes equal to the fidelity of the optimal estimation of thc state ip from N copies 
(Massar and Popescu [1995], Bruss and Macchiavello [1999], Hayashi, Hashimoto 
and Horibe [2004]) 

N+1 

= (4-10) 
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Consequently, in the limit M — > oo the optimal cloning becomes equivalent to the 
optimal State estimation from N copies of ip foUowed by preparation of infinitcly 
many copies of the estimated state. This relationship between optimal universal 
cloning and optimal state estimation can be explored to prové the optimality of the 
cloning transformation (4.4). It foUows from the symmetry, isotropy and linearity 
of universal quantum cloning that the single-qudit outputs must have the form (4.7) 
and that for concatenated universal cloners the shrinking factors multiply. Since 
the concatenation of the optimal N ^ M and M ^ L cloners cannot be better 
than the optimal N ^ L cloner, we get 

7(Ar,L)>7(7V,M)7(M,L). (4.11) 

We take the limit L oo and take into account that the shrinking factor corre- 
sponding to the íidelity (4.10) reads 'y{N,oo) = N/{N + d). In this way, we get 
from the inequality (4.11) an upper bound on 'y{N, M), 

7(iV,M)<2Í^ = ^^, (4.12) 
' ' - 7(M,oo) N + d M ' ^ ' 

which is saturated by the optimal universal cloning transformation (4.4). 



4.1.2. Unitary realization and quantum circuit 

So far the optimal cloning transformation was presented in the form of the rather 
abstract CP map (4.4). It holds that every trace preserving CP map admits a 
unitary realization with the use of an ancilla system. The unitary realization of 
cloning requires 2(M — A^) ancilla qudits; M — N blank copies and M—N additional 
ancillas. For the sake of presentation simplicity we will consider here the N ^ M 
cloning of qubits (Gisin and Massar [1997]). The unitary cloning transformation 
can be expressed in a covariant form, 

M 

[/|ArV>)i„|fí)a„c = ^a.-KM - i)V, ) dones |(M -N- i)V^,i^)anc· (4.13) 

Here \ktp, {N — k)ijj-^) denotes a symmetric state of N qubits with k qubits in state 
IV') and N — k qubits in an orthogonal state IV'"'")) (V'IV'"'") = 0) |-R)anc denotes the 
initial state of the ancilla qubits and 

YM- A^^V M + i\~^/^ 
"^"f-^K n ) [m-n) ■ 

The generalization of the formula (4.13) to qudits with arbitrary d was obtained by 
Fan, Matsumoto and Wadati [2001]. In Sec. 5.1 we shall show that the transfor- 
mation (4.13) naturally arises in stimnlated amplification of light when the qubits 
are represented by the polarization states of single photons. 

Quantum Information theory teaches us that an arbitrary unitary operation 
U can be implemented as a sequence of single-qubit rotations and a two-qubit 
controUed-NOT gates, í7cnot = \j)c\k)t = \j)c\k(Bj)t , where ® denotes addition 
modulo 2, c is the control qubit and t is the target qubit. The quantum network 
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Figure 3: Network for optimal universal 1 M cloning of qubits. The coning 
requires M — 1 blank copies ai, • • • , a^-i and M — 1 ancilla qubits òi, • • • , Òm-i- 
The cloning consists of a preparation of an entangled state |$)afe of blank copies 
and ancillas foUowed by a sequence of 2(M — 1) C-NOT gates between the input 
qubit and the blank copies and ancillas. 



for the optimal universal 1 ^ M cloning of qubits (Buzek, Braunstein, Hillery and 
Bruss [1997], Buzek and Hillery [1998b]) is depicted in Fig. 3. First, the 2(M - 1) 
ancilla qubits ai, • • • , au-i and 61, • • • , 6m-i are prepared in an entangled state 



1 / 2 

* fc=o 

where = M — k and /fe = y/k{M — k), and |M — 1, A:) denotes a symmetric state 
of M — 1 qubits with k qubits in state |1) and M — 1 — k qubits in state |0). The 
state (4.15) can bc gcncratcd by a sequence of singic-qubit rotations and C-NOT 
gates starting from any initial pure state of the ancilla. The cloning itself consists 
of a sequence of 2(M — 1) C-NOT gates where the qubit üq that contains the state 
to be copicd serves as a control qubit and the ancillas are target qubits. This is 
followed by another sequence of 2(M — 1) C-NOT gates where now the qubit üq is 
target and the ancilla qubits are controls. The M clones are stored in the qubits 
dOï • ■ ■ I O'M — 1 while the qubits bj represent the ancillas. 

4.2. OPTIMALITY PROOF FOR 1 ^ M CLONING OF QUBITS 

The optimality of the 1^2 symmetric cloning machine for qubits was first proved 
by Bruss, DiVincenzo, Ekert, Fuchs, Macchiavello and Smolin [1998]. The optimal- 
ity of the cloning transformation (4.4) for arbitrary number of inputs N, outputs 
M and dimension d was proved by Werner [1998] for the global fidelity and later 
by Keyl and Werner [1999] for the single-clone fidelity using powerful techniques 
of group thcory. Hcrc we shall present a simple optimality proof for the class of 
1 ^ M universal symmetric cloning machines for qubits. This proof foUows the 
general concept outlined in Sec. 3.3 where it was shown that the fidelity is upper 
bounded by the maximum eigenvalue of a certain positive semidefinite operator. 
This optimality proof when single-clone fidelity is used as a figure of mèrit is due 
to Gisin and Massar [1997] and it has also been extended to the global fidelity 
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(Fiuràsek [2001b]). The advantagc of this approach is that it can be easily gener- 
alized to asymmetric cloning, as discussed in the next section. 

Consider thc maximization of the single-clone fidelity and let us assume that 
the output Hilbert space of the cloning map S is the symmetric subspace, since 
the desired outputs e Tíf^. Then ah the clones have the same fidehty by 

construction and we can express the operator S that is isomorphic to the CP map 
<S as foUows, 

1 M 

S=Y. 12 Sik,Ji\i)in{j\ \M,k)out{M,l\. (4.16) 

i j=0 k,l=0 

The mean single-clone fidelity can be calculated by averaging over the surface of 
the Poincaré sphere, 

FNtÚ"" = / Tr[(V^ ® V)TrU,(S)]dV, (4.17) 
where Tr^^j denotes tracing over aU output qubits except for the first one, and 

ídil> = ^í í sini?dí?#, 1^) = cos ^|0) + 6*"^ singli). (4.18) 

J4, 47r Jo Jo 2 2 

After the tracing and integration, we find that the mean single-clone fidelity is a 
linear function of S, F^^^'" = Tr[SRsc]) where the positive semidefinite operator 
Rsc reads 

1 ^ 

^sc = ^5][(2M-/c)|0)(0| + (M + /c)|l)(l|]®|M,fc)(M,fc| 

fc=0 

M-l 

+ ^ E V{M-k){k + l)\l){0\ ®\M,k + 1)(M, k\ + h.c. 

fc=0 

According to Eq. (3.20) the fidelity F^^^'^ is upper bounded by the maximum 

eigenvalue of Rsc-, -Pjv'^m — 2rsc,max· The matrix Rgc has a block diagonal 
structure and it is easy to show that all thc cigenstates of Rsc have the form 
a|0)|M, /c) -|-/3|1)|M, A:-|- 1). The calculation of the eigenvalues of Rgc thus reduces 
to finding roots of quadratic polynomials and one finds that Rsc has only three 
diffcrcnt eigenvalues, ri = (2M -|- 1)/(6M), r2 = 1/3, and r^, = 1/6. This provides 
an upper bound F^_^^'" < (2M -|- 1)/(3M) which is saturated by the cloning 
machine (4.4). This proves that the machine (4.4) is optimal. 

Similar chain of arguments can be used to demonstrate the optimality of the 
machine (4.4) when the global fidelity is the figure of mèrit. The mean global 
fidelity can be written as F^_^^ = Tr[SRG], where 

Rg= f i^l^i^Zfdi^. (4.19) 
Jtfi 

With the help of Schur's lemma this integral can be easily evaluated and one obtains 
= M^zs^M+i) where nj^, ^ is a projector onto the symmetric subspace of 
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M + 1 qiibits and Ti denotes partial transposition with respect to the first qubit. 
Again, the matrix Rq is block diagonal and its eigenvalues can be easily determined 
analytically. One finds that r^max = 1/(M+1) which implies F'^'^^f < 2/(M+l) 
and this bound is achieved by the cloner (4.4). 

4.3. UNIVERSAL ASYMMETRIC QUANTUM CLONING 

The quantum cloning machines serve as universal distributors of quantum informa- 

tion among scvcral parties. The symmetric cloner divides the information eqnally 
among all M copies but it is also possible to distribute the information unequally. 
A lot of attention has been devoted to the universal asymmetric 1 — > 2 cloning ma- 
chines for qubits (Ccrf [1998], Buzck, Hillcry and Bednik [1998], Niu and Griffiths 
[1998], Cerf [1999], Cerf [2000a]) and qudits (Cerf [2000b], Cerf, Bourennane, Karls- 
son and Gisin [2002]), which produce two clones A and B with different fidelities 
Fa and Fb ■ The optimal asymmetric cloner can be dcfincd as a machinc that for a 
given fixed fidelity of the first clone Fa maximizes the fidelity of the second clone 
Fb- Such machines can íind applications, e.g., in eavesdropping on quantum key 
distribution protocols, whcrc thcy allow onc to invcstigate the tradc-off bctwccn 
the information gained by the eavesdropper and the disturbance observed at the 
receiver's station. 

In terms of the cloning CP map S, the mean fidelities of the two clones can be 
expressed as Fa = Tr[SRA] and Fb = Tr[SRB], where the positive semidefinite 
operators Rj are given by 

Ra= ^Tn ®^a®Ib dij, Rb= ^Tn ®Ia®^b dtp. (4.20) 
Jjp Jijj 

The optimal asymmetric cloning machine should max;imize a convex mixture of the 
mean fidelities Fa and Fb (Fiuràsek [2003], Lamoureux, Navez, Fiuràsek and Cerf 
[2004], Fiuràsek, Filip and Cerf [2005]), 

F = pFa + {1- v)Fb = Tr[SR], (4.21) 

where R = pRa + (1 p)Rb and p is a parameter that controls the asymmetry 
of the cloner. The maximization of F for a given value of p can be equivalently 
rephrased as a maximization of Fb for a fixed value of Fa- After some àlgebra, we 
find 

R = rf(rf^[^in^AB + dp*+,A ®Ib + -P)$+ b ® ^a]- (4.22) 

The maximum eigenvalue of R is d-fold degenerate and the corresponding eigen- 
vector is 

jw; k)=a \^+)AR\k)B + \^+)BR\k)A, (4.23) 

where the coeflicients a, /3 > are some functions of d and p. By proper ly normal- 
izing the eigenstates (4.23) we get 

a2 + /?^ + ^ = l. (4.24) 

The operator S isomorphic to the optimal cloning CP map S is proportional to 
the projector onto the subspace spanned by the eigenstates (4.23). The unitary 
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rcalization of this map requires a single ancilla qudit C and can be written in a 
covariant way, 

IV) ^ a |V)a|$+)bc + /? |V)b|$+)ac· (4.25) 
From this expression we can evaluate the fidelities of the two clones, 

Fa = 1-^P\ FB = l-^a'. (4.26) 

Note that the parameters and 0^ are the so-called depolarizing fractions as 
discusscd in Cerf [1998], Cerf [2000a], Cerf [2000b]. The one-parametric class of the 
optimal universal asymmetric 1^2 cloning machines is characterized by the Eqs. 
(4.25) and (4.26) together with the normalization condition (4.24). 

4.4. UNIVERSAL NOT GATE 

The process of optimal quantum cloning is closely connected to another impossible 
operation in quantum mechanics, the so-called universal NOT gate for qubits. The 
hypothetical universal NOT gate would perfectly reverse any spin-^ state. This 
device should thus produce from the input qubit |V') an orthogonal state |V'-l). 
However, this is impossible, because the transformation \tp) — » \ip±) is anti-unitary. 
More generally one can consider an extended scenario where A'' copies of the state 
are available and the task is to prepare a single copy of the flipped spin IV'.l). 
The best approximation to this forbidden operation was found in Cisin and Popescu 
[1999], Buzek, Hillery and Wcrncr [1999], Buzek, Hillery and Werner [2000]. 

The optimal universal NOT gate Sunot can be made covariant by twirling so 
that the fidelity F = (V'-l|'5unot('!/'®^)I^-l) does not depend on j^) and can be 
written as .Funot = Tr[SuNOTRuNOT], where 

iïUNOT = / [i'^^'f ® ^±di^ = T7^(C/®^ ® /)n+,A,+i([/t®^ ® J) (4.27) 
J·4> iV + z 

see Fiuràsek [2001b]. The unitary operation U provides the link between the states 
IV'*) and IV-i), IV--*) = U\tpA_), U\0) = -|1), Í7|l) = |0). The fidelity of the U- 
NOT gate is bounded by the maximum eigenvalue of íÏunot- Since this operator 
is proportional to a projector, we immediately find ruNOT.max = l/(-^ + 2). The 
dimension of the input Hilbert space Hf^ is d = N +1 and we obtain, 

ÍUNOT = ]\[ ^2 ' (4.28) 

Remarkably, this fidelity coincidcs with the optimal fidelity of the estimation of the 
state from N copies. If we possess an estimate of \tp) then we can also produce an 
estimate of [í/'.l) with the same fidelity, simply by flipping the estimated spin. This 
implies that the optimal U-NOT gate can be realizcd by pcrforming the optimal 
estimation of the state \ip) foUowed by the preparation of the flipped estimated 
state. In this way we can generate arbitrary many approximate copies of |V'-l)j 
all with the same fidelity (4.28). Remarkably, the optimal cloning transformation 
(4.13) simultaneously also implements the optimal approximate U-NOT gate. This 
machine produces M — N approximate anti-clones, which are stored in the ancilla, 
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Figure 4: Optimal universal cloning of photons by means of stimulated parametric 
down-conversion. BS auxiliary beam splitter, Ci and C2 - nonlinear crystals, WP - 
wave plates for input-state preparation, M ~ mirror, PD ~ triggering single-photon 
detector for conditional preparation of iV-plioton Fock state. 



see Buzek, Hillery and Werner [1999], De Martini, Buzek, Sciarrino, and Sias, 
[2002]. 

It is possible to generalize the concept of U-NOT gate to qudits, by noting that the 
state is unitarily equivalent to the state The complex conjugation is well 
defined for any dimension d, and one can look for the transformation that optimally 
approximates the (generalizcd) transposition map ^ ip* = ip'^ . Using similar 
reasoning as before, one can prové that the maximal fidelity of the approximate 
transposition is equal to the fidelity (4.10) of the optimal state estimation from N 
copies, see Fiuràsek [2004]. 



§ 5. Universal cloning of photons 
5.1. AMPLIFICATION OF LIGHT 

In quantum òptics, single photons are very often used as carriers of quantum Infor- 
mation. Photons represent ideal flying qubits; they can be transmitted over long 
distances via low-loss optical fibers and their interaction with the environment is 
very weak so they do not suffer from a significant decoherence. Quantum bits can 
be encoded into single photons in various ways. One natural option is to exploit 
the polarization degrees of freedom and to represent a qubit as a superposition of 
vertically {\V)) and horizontally {\H)) polarized photon, \ip) ~ a\H) + /3\V). An- 
other possibility is to use the so-called time-bin encoding where the photon can be 
located in one of d different time slots (Marcikic, de Riedmatten, Tittel, Scarani, 
Zbinden, and Gisin [2002], de Riedmatten, Marcikic, Tittel, Zbinden, CoUins and 
Gisin[2004]). Such encoding has the advantage that it is not restricted to qubits and 
the photon can thus represent a d-dimensional system with arbitrary d (de Ried- 
matten, Marcikic, Scarani, Tittel, Zbinden and Gisin [2004]). Arbitrary time-bin 
qudits can be prepared using unbalanced Mach-Zehnder interferometer. Time-bin 
encoding has been advantageously used for long-distance quantum key distribution. 

The cloning of the quantum states of single photons requires that the number 
of output photons be higher than the number of input photons. This simple fact 
immediately leads to the insight that the optimal copying of photons can be per- 
formed by means of amplification of light (De Martini, Mussi and Bovino [2000], 
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Simón, Wcihs and Zcilingcr [2000a], Simón, Weihs and Zeilinger [2000b]). This is 
very natural bccausc thc goal of quantum cloning is to "amplify" the quantum in- 
formation carried by thc photons. There are several physical mechanisms that can 
bc uscd for cloning sucli as parametric down-conversion or amplification of light 
in atòmic media. In all cases, the cloning is achicvcd due to the process of stimu- 
lated emission, which means that the medium emits preferably photons in the same 
quantum state as that of thc input photons injected into the medium. 

Most of thc quantum cloning experiments based on the stimulated amplifica- 
tion of light wcre carried out using the process of stimulated parametric down- 
conversion. Consider a nonlinear crystal with sccond ordcr nonlincarity x^"^^ ■ In 
such crystal, a single "blue" pump photon with frcquency u)p can be converted 
into two "red" photons with freqüències and W/ such that + = Wp, 
which expresses the energy conservation. The two downconverted photons are re- 
ferred to as signal (S) and idlcr (I), respectively, for historical reasons. An efScient 
down-conversion requires thc conservation of momentum, which translates into the 
phase-matching condition k^ + k/ = kp, where kj stands for thc wavc- vector of the 
j-th photon and |m — njujj/c, where rij is the refraction indox at frequency ujj. 
Efhcient phase matching in thc nonlinear crystal can be achieved by exploring the 
birefringence and using different polarizations for the pump, signal and idler beams. 
We can distinguish two different kinds of phase matching. In Type I matching, the 
pump beam is, say, vcrtically polarized, and both signal and idler are horizontally 
polarized. On the other hand, in the Type II matching, the signal and idler photons 
are orthogonally polarized. Besides their polarization states, the signal and idler 
beams can also be distinguished spatially. So in non-degenerate Type-II downcon- 
version we deal with modes Ah and Ay for the signal beam and Bh and By for 
the idler beam. It is possible to arrange the configuration of the pump beam and 
nonlinear crystal and to select only certain directions in the output beams such 
that the effective Hamiltonian which describes this process reads 

H = ÍK{alb^H - a'fjblr) + h.c. (5.1) 

This Hamiltonian is obtained in the limit of strong coherent pumping and the 
coupling constant k is proportional to the pump beam amplitude ap and to the 
second-order nonlinearity x^^^ , while is the creation operator for the j-th mode. 
An essential feature of thc Hamiltonian (5.1) is that it is invariant with respect 
to the simultaneous identical transformation of the polarization basis of signal and 
idler photons. Mathematically, we have (U ®U)H{W (g) W) = H, where UavU^ = 
uvvfív + uvHdH, UanU^ = uhvo-v + uhh<JH: thc matrix Uij is unitary and 
identical transformation rules hold also for by and bn- This covariance property 
guarantees that the cloning process is universal and the cloning fidelity is the same 
for all input states. It thcrcforc sufficcs to consider only one particular input state. 

A sketch of the cloning of polarization states of photons via stimulated downcon- 
version is shown in Fig. 4. The signal mode is initially prepared in A''-photon state 
\tp)a^ ■ This can be in practice achieved c.g. by means of spontaneous parametric 
downconversion in crystal Ci and conditioning on observing N photons in the out- 
put idler mode with photodetector PD. After the passage through the crystal C2, 
M — N photon pairs can be generated with certain probability. If this happens, 
then M clones are present in the mode A while the mode B contains M — N anti- 
clones. Note that the cloning is only probabilistic and we cannot predict a-priori 
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Figure 5: Experimental setup for optimal universal cloning by means of stimu- 
lated parametric down-conversion. The input single-photon state to be cloned is 
obtaincd from a weak coherent làser beam. (After Lamas-Linares, Simón, Howell 
and Bouwmeester [2002].) 



the number of clones that will be generated. The particular N ^ M cloning event 
can be selected only a-posteriori by accepting only the events when M photons 
were detected in mode A. 

Let us start with a simple example of 1 — > 2 cloning which will illustrate all the 
main features. In this case, the input state is given by a single photon in mode A and 
a vacuum in mode B. As already explained, without loss of generality we can assume 
that the photon is vertically polarized and we have jT/iin) = |l)av|0)ajf |0)b^|0)b„. 
The generation of the second clone requires that a single photon pair is emitted 
in the nonlinear crystal. In the first order perturbation theory, the output state is 
given by 

H\l)av\0)a„\0)b.\0)b„ CX V2\2)a,\0)a,\0)bjl)b„ - 1 1) 1 1) a„ 1 1) 6 J 0) fc„ . (5.2) 

Notice the prefactor V2 which arises because the emission of the second vertically 
polarized photon in mode A is stimulated by the presence of a vertically polarized 
photon in this spatial mode. This cloning is optimal since it yields the maximum 
fidelity. We can immediately see that the global fidelity is 2/3. To determine the 
single-clone fidelity we note that with probability 2/3 both photons in mode A are 
vertically polarized and with probability 1/3 only one photon is vertically polarized. 
So the probability that one randomly chosen photon in spatial mode A is vertically 
polarized is|xl + -|xi = |, which is the maximal single-clone fidelity for 1 — > 2 
cloning of qubits. 

Several experiments on cloning via parametric down-conversion have been re- 
ported, see De Martini, Mussi and Bovino [2000], Lamas-Linares, Simón, Howell 
and Bouwmeester [2002], Pelliccia, Schettini, Sciarrino, Sias and De Martini [2003], 
Sias, Sciarrino and De Martini [2003], De Martini, Pelliccia and Sciarrino, [2004]. 
The experimental setup used by Lamas-Linares, Simón, Howell and Bouwmeester 
[2002] is shown in Fig. 5. A nonlinear BBO crystal is pumped by a second harmònic 
of a Ti:sapphire làser which emits 120 fs long pulses. A tiny part of the coherent 
master làser beam was split on the first beam splitter BS and used as a seed for the 
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Figure 6: Observed coincidence rates in dependence on the position of the movable 
mirror of two clones in identical polarization states (A,B,C) and in in orthogonal 
polarization (D,E,F) for three different input polarizations. (After Lamas-Linares, 
Simón, Howell and Bouwmeester [2002].) 



down-conversion. With probability p <^ 1, the beam contained exactly one photon. 

This beam was fed to the BBO crystal and the output is analyzed using a sequence 
of wave-plates, polarizing beam splitters (PBS) and single-photon photo-detectors. 
The probabiUty of a pair generation in the crystal P2 <^ P which guarantees that 
the dominant event leading to two photons in mode a and one photon in mode b 
is when a single photon was in the weak coherent beam a and a single pair was 
emitted in the crystal. The conditioning on observing a click of the trigger detector 
Dl is important since it eliminates the events when there were two photons in mode 
a and no pair was generated in the crystal. 

In the experiment, one measures the number of coincidence clicks of the pho- 
todetectors D2 and D3 as a function of the time delay between the input photon 
beam in mode a and the pump beam. If those two beams do not overlap in the 
BBO crystal, then there is no stimulated down-conversion and the polarization of 
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the second photon cmitted in mode a is fuUy random. If the two bcams overlap, 
then stimulated amplification sets on and the second photon is emitted preferably 
with the same polarization as the input photon. Optimal cloning is achieved when 
the overlap is perfect. A detector setting with PBS was used to measure the num- 
ber of orthogonally polarized photon pairs iV(l, 1). To detect the number of pairs 
with the same polarization iV(2, 0), PBS was replaced by a polarizer followed by an 
ordinary beam sphttcr. The observed coincidence rates as a function of the time 
delay arc shown in Fig. 6 for three difFerent polarizations. We can sec that N{í, 1) 
docs not depend on the delay as expected, whilc iV(2, 0) decreascs with increasing 
dclay. The avcrage experimental cloning fidclity dctcrmincd from thcse data reads 
F « 0.81 which is very close to the theoretical maximum 5/6 « 0.833. 

An improved experimental setup involving double passage of the pump beam 
through the nonlinear crystal has been developed (Pelliccia, Schettini, Sciarrino, 
Sias and De Martini [2003], De Martini, Pelliccia and Sciarrino, [2004]), sec Fig. 7. 
In this setup, the photon to be cloncd is gcncrated during the first passage of the 
pump pulse through the crystal. Sincc the signal and idler beams are entangled, 
projecting idler beam onto statc j^) prepares the signal in state {ip^). The click of 
the trigger detector Dt heralds the preparation of a single photon in mode labeled 
— fci in Fig. 7. This photon is then cloned by sending it again through the nonlinear 
BBO crystal. The delay between the pump and signal is controUed by moving the 
mirror Mp. In this experiment, the states of both clones and the anti-clone were 
analyzed simultaneously and it was demonstrated that this device accomplishes 
jointly the optimal 1^2 cloning and also the optimal universal-NOT gate for 
qubits. The attained fidelities were Fclon = 0.81 and Funot = 0.62. 

The stimulated down conversion can be used to probabilistically implement any 
N ^ M cloning of qubits (Simón, Weihs and Zeilinger [2000a]) and even qudits 
(Kempe, Simón and Weihs [2000], Fan, Weihs, Matsumoto and Imai [2002]). The 
unitary transformation induced by the Hamiltonian H can written in a factorized 
form as foUows, 



e 



where A = tanh(ííí), t is an effective interaction time, and ntot = üyav + a\jaH + 
bybv + b^H^H is the total number of photons in spatial modes a and b. Since 
the Hamiltonian H is covariant it is enough to consider the input state \ipin) = 
\^)av\0}aH\0)bv\0)bH- With the help of the factorization (5.3) we find that the 
corresponding output state reads 

oo 

e-'"'\An) = (1 - E A^-^'I^m), (5.4) 



M=N 



where 



M-N I . _ , s 

l*^)=E(-l)y( N )\M-k)av\k)aH\k)bv\M-N-k)m- (5.5) 

We can see that the output state (5.4) is a weighted superposition of states \^m) 
with difïerent numbers of clones M. The state \^m) and hence the fidelity of 




Figure 7: Simultaneous contextual realization of the optimal cloning machine and 
the universal NOT gate. Q - quarter wave-plates, WP - half-wave plates, PBS - 
polarizing beam splitters, BS - beam splitters, BBO - nonlinear crystal, D ~ single- 
photon detectors, M - movable mirror. (After De Martini, Pelliccia and Sciarrino, 
[2004].) 



N M cloning is independent of the coupling strcngth A and only the probability 
of generating exactly M clones depends on A. One can also inimediately see that 
the state (5.5) coincides (up to an irrelevant overall normalization factor) with the 
outcome of the optimal cloning transformation (4.13), hence the universal cloning 
via parametric down-conversion is optimal. 

We now extend the concept of cloning via amplification to qudits represented 
by a single photon in d different spatial modes or time-bins. The use of time-bin 
encoding seems to be particularly advantageous since only a single nonlinear Type-I 
matched crystal is required, and the pump beam should consist of a sequence of d 
pulses. We associate creation operators aj and bj with j-th time-bin of signal and 
idler beams, respectively. The Hamiltonian governing the evolution of this system 
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can be expressed as 

d 

Ha = ÍKj2{a]b]-ajbj). (5.6) 

j=i 

This Hamiltonian is invariant with rcspcct to simultaneous unitary transforniations 
of the signal and idler modes, {U ® U*)H{W (g> U^) = H, where U G SU{d). 
This covariance property guarantccs that thc cloning is universal and thc cloning 
fidehty does not dcpcnd on thc input statc so it sufhccs to considcr thc input state 
IV'in.d) = \N)ai |0)a2 ' ' ' |0) a J 0) J 0) ' ' ' 10)6^. Thc unitary operation cxp(-iiíí) 
can be again factorized siniilarly as in Eq. (5.3) and wc gct e-'"''^\ip;nA) = (1 - 
^2)iv/2+d^~^^ A^-^lí-M.d) where the state containing M clones reads 



I^M.d) = XI y iv"^^) "^^^"^ \'^·i)o.^ ■ ■ ■ \i^d)ai \mi)b^ \m2)b2 ■ ■ ■ \md)bi ■ 

(5.7) 

In this formula, indicates summation ovcr all vectors m = (toi, . . . , m^) sat- 
isfying Xlj=i ^3 = M — N. The optimality of this cloning transformation can 
be proved by explícit evaluation of the fidelity. It can be shown that there are 
(^~^d^2~^~'") different terms in Eq. (5.7) with N + m photons in mode ai, each 
with weight {^^) ■ The average single-clone fidelity can be thus expressed as. 



1 ^í-^^ fN + m\fM-N + d-2-m\N + m 



E 

m=0 

where the normalization factor is given by 

M-N 



^ I N + m\ fM - N + d-2-m\ fM + d-ï\ 
E M d-2 = Ar^.-i • (5-9) 



m=0 



The summation in Eq. (5.8) can be performed with the help of the identity given 
in Eq. (5.9) and onc rccovcrs thc optimal fidelity (4.9). 

Instead of parametric down-conversion it is also possible to amplify the light by 
sending it through an inverted atòmic medium (Simón, Weihs and Zeilinger [2000a] , 
Kcmpc, Simón and Weihs [2000], Fan, Wcihs, Matsumoto and Imai [2002]). Thc 
àtoms should possess d diS'erent ground states \gj) and an excited state |e). We 
assume that each atòmic transition |e) — > \gj) is strongly coupled to a single optical 
mode aj and thc qudits are rcprcscntcd by single photons in thosc d modes. Thc 
universality of the cloning requires that the coupling strength k must be the same 
for all d transitions |e) — > [pj). In the interaction picture and in the rotating 
wave approximation, the interaction of light with àtoms is governed by the Jaynes- 
Cummings Hamiltonian, 

L d 

Hjc = KX5I4l5í·ít)(efc| +h.c., (5.10) 

fc=ii=i 

where L is thc numbcr of àtoms and \gjk) stands for the ground state \gj) of k-th 
atom. This Hamiltonian satisfies the covariance property U (8) U*HW (8) U'^ = H, 
where Ua'jW = ujkal, U\ek) = \ek), and U\gjk) = Ylf=iUji\gik). 
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Suppose that all L = M — N àtoms are initially prcparcd in the excited state 
and that all N input photons are in mode Ai. The joint atoms-photons input state 
reads IV'íii.la) = |-^)ai |0)a2 • • • |0)ad|ei) • • • Icl). If each atom emits a photon during 
the passage of the light through the àtoms then M clones are generated and all 
àtoms end up in ground states. 

In the regime of weak coupling we can express the output state conditional on 
all àtoms being in some ground state using the L-th order perturbation theory, 

d 

IV'out) CC (X^a]ò].c)^ |V'in,LA), (5.11) 

i=i 

where the operator 6]c is defined as bjC = J2k=i \9jk){ek\- Note that with this nota- 
tion, the Hamiltonian (5.10) becomes similar to the down-conversion Hamiltonian 
(5.6). Since the àtoms arc supposcd to be idcntical, the photons cmitted by them do 
not carry any Information about which atom emitted which photon. Consequently, 
if all àtoms emit photons, then the àtoms relax to symmetric ground state. Sup- 
pose that iTij photons were emitted to mode aj, j = l,...,d. The corresponding 
symmctrizcd ground atòmic state reads 

Iffm) = ^ \giki) . . . \gik,ni)\92k,m + l) ■ ■ ■ \92k,m+m2) ■ ■ ■ \9dkL-ma + l) ■ ■ ■ \9dkL)^ 

7r(k) 

(5.12) 

whcrc X^^çk) denotes summation ovcr all L\ vàlues of the subscripts fc;, í = 1, . . . , L, 
which can be obtained as permutations of {1, . . . , L}. The normalization coefhcient 

CÍ^ = mi!m2!···md!L! (5.13) 

is choscn such that (.gm|.9m) = 1- After some àlgebra, one finds that the output 
state (4.6) can be expresscd as follows, 




+ mi)ai l"Ï2)a2 • ■ • \'md)ai\9m)· 

Since this state is fuUy equivalent to the state (5.7), the cloning is optimal. Although 
this result was obtained within the framework of a perturbation theory, a detailed 
analysis reveals that it holds for any interaction strength. It can be also shown 
that if only M' — N < L àtoms emit photons and the rest of the àtoms remain 
in the excited state, then M' optimal photonic clones are generated (Fan, Weihs, 
Matsumoto and Imai [2002]). 

A proof-of principle experiment on cloning via stimulated emission was reported 
by Fasel, Gisin, Ribordy, Scarani and Zbinden [2002] . In that work, a commercially 
available polarization insensitive Erbium doped fiber amplifier was utilized. The 
amplifier was injectcd with a weak vertically polarized coherent signal with mean 
photon number fiin- After the ampliíication, the output mean numbers fïy and üh 
of vertically and horizontally polarized photons were measured. The íidelity of the 
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amplification process can be simply defincd as F = nv/{ny + nu). Tlic output 
mean intensities depend linearly on the input intensity (Shimoda, Takahasi and 
Townes [1957]), 

ny = Gnin + ^(G-l), njï = ^(G-l). (5.14) 

Here G is the gain of the ampHfier and Q is a factor depending on thc propcrties of 
the ampUfication process. Thc term Güi^ representa the amplified injected input 
signal while (G — 1)/Q represents thc noise arising due to spontaneous emission. 
For quantum- noise Umited amplification, Q = 1. From Eqs. (5.14) we can express 
G in terms of nin, fïout = nv+fiH and Q, and we find G = (Qíï-out + 2)/((5nin + 2). 
On inserting this expression into thc formula for the fidehty, we obtain 

p Q^out^in ~t~ ?^out "i" ^in \K\ 

Qnontfiin + 2nout 

If we formally replace üin with A'' (the number of input copies) and üout with M 
(the number of output clones), then for Q = 1 thc formula (5.15) bccomcs the 
optimal fidelity oí N M cloning of qubits. Experimentally, G = 1.3 and Q = 0.8 
was observed which is quite close to the optimal value Q = 1. For instance, the 
fidelity of 1 — > 2 cloning for Q = 0.8 inferred from (5.15) reads F = 0.821 which is 
only slightly lower than the optimal fidelity f = 5/6 » 0.833. 



5.2. SYMMETRIZATION 

We have seen in Section 4.1 that the optimal universal — > M quantum cloning 

can bc accomplishcd by symmctrizing thc statc of N input copies and M — N 
maximally mixed states. Since photons are bosons, the projection onto symmetric 
subspace can be easily carried out with the use of linear òptics, namely by mixing 
thc Aí photons on an array of M — 1 bcam splittcrs and selecting only the events 
when all photons are coUected in a single spatial mode. 

Let us first illustrate this method on the example of 1 — > 2 cloning of polarization 
states of photons (Ricci, Sciarrino, Sias and Dc Martini [2004], Irvinc, Lamas- 
Linares, de Dood and Bouwmeester [2004], Sciarrino, Sias, Ricci and De Martini 
[2004b]). The setup is schematically illustrated in Fig. 8(a). The photon in mode 
A whosc statc is to be cloncd is combincd on a balanccd bcam splittcr BSi with 
a blank copy photon prepared in a maximally mixed state. Only the cases when 
both photons leave the beam splitter in the left output mode are post-selected and 
the two clones arc spatially scparatcd by an auxiliary balanccd bcam splitter BS2· 
At the heart of cloning via symmetrization is the Hong-Ou-Mandel effect (Hong, 
Ou and Mandel [1987]). If two photons with identical polarization state interfere 
on a balanccd bcam splitter, thcn thcy both cnd up in thc samc spatial mode and 
one does not observe any coincidences of one photon in mode A and one in mode 
B. So for the input there is probability 1/2 of having two photons in the 

left output port and probability 1/2 of splitting thcm in the two output modes A' 
and B'. Altogether, the conditional transformation reads 



I^)a|V')s-^IV')a'|^)b', 



(5.16) 



5. UNIVERSAL CLONING OF PHOTONS 



51 



pump 




(b) 

l¥> 











^ BS'i 










— ? 

BSi 






BSm-1 






H 

; : 


1 ■■■■ \ 





|\|/> |\|/> p=I/2 p=I/2 ' 

N input States M-N blank copies 



M clones 



PD 



Figure 8: Optimal cloning of polarization statcs of photons via projection onto sym- 
metric subspace. (a) The optimal 1-^2 universal cloning based on the interference 
of two photons on a balanced beam splitter BSi. (b) Extension to the optimal 
N ^ M cloning. The A^" input states and M — N blank copies in maximally mixed 
States are combined on an array of M — 1 beam splitters BSj and the clones are 
then separated on another array of M — 1 beam splitters BS'j. 



On the other hand, if the two photons are initially in orthogonal polarization states, 
|V')a|V'_l)_b, then they are distinguishable and do not interfere on BSi. With prob- 
ability 1/4, the photon in mode A is reflected and the photon in mode B is trans- 
mitted and they are both in the left output. Again, there is probability 1/2 that 
the two photons will be divided on a balanced beam splitter BS2. Since the photon 
in state can be either reflected or transmitted on BS2, the final state of photons 
in modes A' and B' is a balanced superposition of these two possibilities, namely 
a symmetric state, 

IV')a|V'í)b - \{\^)a'\Í^±)b' + \Í^±)a'\Í^)b'). (5.17) 

Since the projector onto the symmetric subspace acts as n+j-í/))!?/)) = ^-nd 
n_(_|'0)|'0j^) = + it immediately foUows from Eqs. (5.16) and 

(5.17) that the setup shown in Fig. 8(a) implements with probability 1/4 the pro- 
jection onto the symmetric subspace foUowed by a spatial separation of the two 
photons. 

The maximally mixed polarization state in mode B can be obtained for instance 
by preparing the blank copy photon in state \V) or \H) with probability 1/2 each. 
Another, moro intriguing, option is to send into the port B one part of the maximally 
entangled two-photon singlet state jí'^) = ^(|V')b|V'-l)c — |V'-l)bV')c)· In this case, 
if the symmetrization succeeds, then we obtain in the spatial mode C the optimal 
anti-clone of i.e. a state that has a fidelity 2/3 with 

The optimal 1 — > 2 cloning based on symmetrization has been experimentally 
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demonstratcd by two groups (Ricci, Sciarrino, Sias and De Martini [2004], Irvine, 
Lamas-Linarcs, dc Dood and Bouwmeester [2004]). In both experiments, the input 
photon whose statc was cloned was obtained from a weak coherent beam, and it was 
combined on a balanccd bcani splittcr with one photon from a maximally entangled 
singlet State generated in a nonhnear crystal by means of spontaneous parametric 
down-conversion. The triple coincidence events were selected where there were 
two clones and one anticlone present and the intensity of the weak coherent beam 
was adjusted such that the dominant contribution to the triple coincidence events 
originatcd from the cases when there was a single photon in the coherent beam and 
a singlo entangled photon pair was generated in the nonlinear crystal. The observcd 
mean cloning fidelities in these two experiments were F ~ 0.82 (Ricci, Sciarrino, 
Sias and De Martini [2004]) and F = 0.81 (Irvine, Lamas-Linares, de Dood and 
Bouwmeester [2004]), respectively. The simpler setup depicted in Fig. 8(a) and 
involving only two photons was also implemented experimentally (Sciarrino, Sias. 
Ricci and De Martini [2004a]). A single photon pair was generated in a nonlinear 
crystal. One photon representing the input was prepared in the state [ip) using 
wave plates whilo the other photon was randomly prepared in the state \V) or \H). 
This experiment is much simpler than the previous one, because only two-photon 
coincidence events were observed instead of tree-photon coincidences. This resulted 
in much higher rate of cloning, and also in better visibility and mean cloning fidelity 
F = 0.826 very close to the theoretical maximum F = 0.833. 

An extension of the symmetrization procedure to M photons is illustratcd in 
Fig. 8(b). The photons are combined on an array of M — 1 beam splitters BSj 
and the symmetrization succeeds if all M photons are bunched in the same spatial 
mode (Sciarrino, Sias, Ricci and De Martini [2004b]). To confirm this we can split 
the output signal into M different spatial modes using another array of M — 1 
beam splitters BS^- and postselect only events when each of M photodetectors PD 
registers one photon. 

We now demonstrate that the array of beam splitters accomplishes the desired 
projection onto the symmetric subspace. The symmetric two-mode i-photon states 
|L, k) with L — k photons polarized vertically and k photons polarizcd horizontally 
form a basis in the symmetric space Tí^^ . We prové our claim by induction. 
Consider the Lth beam splitter BSi, in the scheme in Fig. 8(b). The state impinging 
from the left is a symmetric L-photon state while a single photon impinges on BSl 
from the bottom. The beam splitter BSl need not be balanced but its transmittance 
t and reflectance r should be independent of the polarization. In the Heisenberg 
picture, the mixing of the modes on the beam splitter is described by linear input 
output canonical transformations of the creation operators, 

and similar formulas hold also for horizontal polarization. The state transformation 
on a beam splitter can be most easily dctermined by expressing all states in terms 
of the creation operators acting on the vacuum, 

1^' ^> = /i.ul Ml <in<in Vac), |V) = b^.ijvac), |H) = bjj .Jvac). 

(5.19) 
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From Eqs. (5.18) we express the "in" operators as linear combinations of thc "out" 
operators and substitute into the formulas (5.19). Using this technique it is easy to 
show that if all L + 1 photons bunch in the right output mode then the following 
conditional transformation takes place: 



\L,k)\V) tÍrLVL+í-k\L+l,k), 

\L,k)\H) ^ tírLVk + ï\L + l,k + l). (5.20) 

Considcr now thc projcction of thc statcs |L, k)\V) and \L, k)\H) onto the symmetric 
subspace of L + 1 photonic qubits. One finds that 



n+,i+i|L,fc)|y) = yí±l^|L + i,fc), 



n+,i+i|X, k)\H) = ^ ^±^\L + l,k + 1). (5.21) 

The transformations (5.20) and (5.21) are equivalent up to state-independent pref- 

actor \/L + \ tj^r^ which proves that the array of M — 1 beam splitters BS in Fig. 
8(b) projects thc input states onto the symmetric subspace of M qubits. The prob- 
ability of success of the projection can be determined by comparing the coefRcients 
in Eqs. (5.20) and (5.21) and we find 

M- 1 

P = T/(l-T,), (5.22) 

where Tj = í| and Ps = Tr[n^pi„] is the overlap of the input M-photon state pin 
with projector onto the symmetric subspace. Thc optimal transmittance Tj of the 
jth beam splitter leading to maximal P can be obtained by maximizing Tj{l — Tj) 
which yields Tj^opt = + !)• Note that Tj.opt does not depend on the input 
A''-photon state. On inserting the optimal Tj into Eq. (5.22) we get 

M! 

Po,t = Psj^. (5.23) 

Recently, thc optimal universal 1^3 and 2^3 cloning of polarization statcs 
of photons via symmetrization has been demonstrated experimentally in MasuUo, 
Ricci and De Martini [2004]. The three photons used in the experiment consisted 
of a pair of photons gcncratcd in the process of spontaneous parametric downcon- 
version and a single photon in very weak coherent beam. These three photons were 
combined on two beam splitters and only the events where all photons bunched 
in a single spatial mode were chosen by post sclcction. Wave plates, a polarizing 
beam splitter, an array of beam splitters and photodetectors were employed to an- 
alyze the clones. The experimentally observed fidelity of the 1 — > 3 cloning was 
^1^3 — 0-758 which is very close to the theoretical maximum F^'^g = 7/9 ~ 0.778. 
The observed fidelity of 2 ^ 3 cloning, ^"3^3 = 0.894, was also close to the optimum 
value F^ts = 11/12 « 0.917. 

The symmetrization on a beam splitter can be naturally extcndcd to qudits. 
Symmetrization of two photonic qudits represented by a state of a photon in d 
difïerent spatial modes would require an array of d balanced beam splitters, each 
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mixing the j-th mode of thc first and second qudits. It may be more advantageous 
to work with time-bin qudits, where the symmetrization would requirc only one 
balanced beam splitter where the two photons would interferc. Similarly as before, 
only the evcnts when the two photons bunch and leave the beam splitter in the 
same spatial port have to be post-selected. 



5.3. UNIVERSAL ASYMMETRIC CLONING OF PHOTONS 

So far, we havc prcscntcd various optical implcmcntations of symmctric cloning 
machines. In this section we will consider optimal 1^2 asymmetric cloning of 
qubits. We will describe two methods, both based on the interference of photons 
on unbalanccd beam splittcrs. The first approach, introduccd in Section 3.5, is to 
start from the output of the optimal symmetric cloner and convert it into an output 
of the optimal asymmetric cloner, which is given by 

I*) = , j ^ , [|^)a|V')b|^±>c -p|V')a|^±)b|^)c - {1-pM±)a\Í^)bWc] 
V2 - 2p + 2p^ 

(5.24) 

Here p G [0, 1] is an asymnictry parameter and the fidelities of the clones in qubits 
A and B read 



The symmetric cloner is recovered when p =1/2 and we have 



= [2|^)^|^)b|^^)c - |V)a|V±)b|^)c - |^±)a|V)s|V')c] . (5.26) 

Suppose first that the second clone (qubit B) and the anti-clone (qubit C) are 
projected on the singlet state l^*"). We obtain 

/^®nsc;|í')sym = ^|^)A|*-)BC, (5.27) 

where = The original input state IV») is perfectly recovered in the 

qubit A. The projection on singlet forms a part of the Bell measurement, i.e. a 
measurement in the basis of four maximally entangled Bell states. There is an 
interesting analogy between Eq. (5.27) and the process of quantum teleportation 
(Bennett, Brassard, Crepeau, Jozsa, Peres and Wootters [1993], Bouwmeester, Pan, 
Mattle, Eibl, Weinfurter and Zeilinger [1997], Boschi, Branca, De Martini, Hardy 
and Popescu [1998], Marcikic, de Riedmatten, Tittel, Zbinden and Gisin [2003]). 
In quantum teleportation, two parties, Alice and Bob, share a maximally entangled 
singlet state. Alice wants to send to Bob an unknown state ji/;) using only classical 
communication and shared entanglement as a resource. To achieve this, Alice 
carries out a Bell measurement on the state and her part of the shared singlet 
state. She communicates the measurement outcome (two classical bits) to Bob, 
who accordingly applies one of four different unitary transformations to his part of 
the singlet. In this way, the state is teleported from Alice to Bob. 

The cloning can be deterministically reversed by performing Bell measurement 
on one of the clones and the anti-clone and applying an appropriate correcting 
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Figure 9: Optimal universal asymmetric cloning from symmetric cloning. The 
input photon to be cloned is obtained from a weak coherent beam. BS - auxihary 
beam sphtter, SHG - second harmònic generation, C - nonhnear crystal, 50 : 50 - 
balanced beam spHtter, T : R - unbalanced beam sphtter, WP - wave plates for 
input-state preparation, M - mirrors. 



unitary to the first clone (Bruss, Calsamiglia and Lütkenhaus [2001]). In the Bell 
measurement, the singlet is detected with probability 3/4 while each of the triplet 
Bell states is detected with probability 1/12, independently of the input state. 

This full reversal of cloning can be generalized to a partial reversal which converts 
the symmetric cloner to asymmetric one (Filip [2004a]). The idea is to apply to 
qubits B and C a filter + all^p, whcrc a e [0,1] controls the asymmetry. 

If a = 0, then we get projection on singlet and full reversal, while for a = 1 the 
two qubits are multiplied by identity and nothing happens. Let us now consider 
arbitrary a. The state after filtering, 

lí-proj) ^Ia® {^BC + «n+c.)|í'sym), (5.28) 

can be after normalization expressed as foUows, 

jM/proj) = , ^ \ [(3 + amA\^)B\^i^c - (3 - amA\^i^B\i>)c 
V6(a^ + 3) 

-2a\i^i)A\^)Bmc] (5.29) 

We can immediately sec that this state coincides with the outcome of the optimal 
asymmetric cloner (5.24) and p — {"ò — a)/(3 + a). 

For optical polarization qubits, the filtration (5.28) can be implemented by letting 
the two photons interfere on an unbalanced beam splitter and post-selecting only 
the events when a single photon is detected in each output port. There are two 
ways how the photons can exit the beam splitter in different spatial modes. Either 
both photons are reflected or they are both transmitted. The unitarity dictates 
that these two alternatives acquire a mutual phase shift tt. If the two photons are 
in the same state j·í/'), then these two alternatives interfere destructively, while if 
the photons are in orthogonal polarization states then there is no interference. The 
resulting conditional transformation reads 

\^Í>)bc^{R-T)\í:^)bc, IV'V'i) ^iï|^^±)-T|V'iV)· (5.30) 

It foUows that the unbalanced beam splitter applies the filter 11^ + «11+ with 
a = i? — T. A schematic setup of the proposed asymmetric cloning experiment 
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Figure 10: Partial quantum state teleportation as optimal asymmetric cloning. The 
input photon to be cloned is obtained from a weak coherent beam. BS auxihary 
beam sphtter, SHG - second harmònic generation, Ci - nonhnear crystal, T : R 
— unbalanced beam splitter, WP - wave plates for input-state preparation, M - 
mirrors. 



is shown in Fig. 9. Optimal symmetric cloning is accomplished by stimulated 
parametric downconversion as discussed in detail in Sec. 5.1. At the output, the 
two clones are separated on an auxiliary balanced beam splitter and one of the 
clones is combined with the anti-clone on an unbalanced beam splitter. Successful 
asymmetric cloning is heralded by a coincident observation of a single photon in 
each mode A, B, and C. 

The second scheme for optimal asymmetric cloning (Filip [2004b]) resembles very 
closely the scheme for teleportation of polarization states of photons, see Fig. 10. 
The only difference is that the balanced beam splitter which is used in teleporta- 
tion to perform a Bell analysis is replaced with an unbalanced beam splitter that 
conditionally applies the filter n_ + «11+ . The cloning succeeds if a single photon 
is detected in each mode A, B, and C. The initial state in the scheme shown in Fig. 
10 is \iP)b\'^~)ac and after the interference on a beam splitter and postselection 
we get 

lí-proj) =Ia® iaU+c + UgcMB\'í>~)AC (5.31) 

After some àlgebra we arrive at 

lí-proj) (X \ [il + aMAWB\i^±)c~a~amA\^^)B\^)c 

-2a\^A_)AH)B\Í')c]- (5.32) 

This is again the output state of the optimal asymmetric cloning machine with 
p = (1 — a)/(l + a) so the asymmetric cloning can be implemented by means of a 
partial teleportation. An interesting feature of this scheme is that one of the clones 
is teleported from Alice to Bob so we can speak about cloning at a distance. 

The universal asymmetric cloning of polarization states of single photons has 
been demonstrated experimentally by Zhao, Zhang, Zhou, Chen, Lu, Karlsson and 
Pan [2005] foUowing the scheme illustrated in Fig. 10. In that experiment, a 
Mach-Zehnder interferometer acted as an effective unbalanced beam splitter whose 
transmittance could be controUed by changing the relative path difference between 
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the two arms of the interferometcr. In this way it was possible to demonstrate the 
whole class of the asymmetric 1^2 cloning machines. 

5.4. CLONING OF ORTHOGONALLY-POLARIZED PHOTONS 

It was shown in Scction 4.1 that the optimal universal quantum cloning and optimal 
quantum state estimation are closely related and that in the limit of infinite number 
of clones the fidelity of cloning is equal to the fidelity of optimal state estimation. In 
this context, a very intcrcsting and snrprising obscrvation was madc by Gisin and 
Popescu [1999] , who found that the state of a single qubit can be better estimated 
from the state than from the state |V')IV')· Picturing the qubits as spin-i 

partides, we can say that the Information about the dircction is bcttcr cncodcd 
in two anti-parallel spins than in two parallel ones. The fidelity of the estimation 
of IV') from a single copy of the two-qubit state |V')|V'-l) reads (Gisin and Popescu 
[1999], Massar [2000]), 

= ^1 + ^) « 0.789, (5.33) 

which is slightly higher than the fidelity of optimal estimation from |V')IV')) = 3/4. 
Motivated by this observation wc may expect that this advantagc of two anti- 
parallel spins over two parallel ones extends also to the cloning. This is indeed the 
case, provided that the number of clones is large enough. We shall now describe the 
optimal universal cloning transformation which produccs M approximate clones of 
the state IV») from a single replica of |V')|V'-l) and maximizes the single-clone fidelity. 

Making the natural assumption that the output Hilbert space is the symmetric 
subspace of M qubits, the optimal cloning CP map S can be dctcrmined analytically 
for any M (Fiuràsek, Iblisdir, Massar and Cerf [2002]). The mean single-clone 
fidelity can be expressed as F = Tr[SR], where R has a rather complicated form 
and can be found in (Fiuràsek, Iblisdir, Massar and Ccrf [2003]). In contrast to 
the universal cloning with input state \ip)^^ , the maximum fidelity cannot be 
determined from the maximum eigenvalue of R, and one has to solve the extremal 
Eqs. (3.18) and prové the optimality by checking that the inequality (3.19) is 
satisfied. 

Since the input state of the cloner can be obtained as an orbit of the group 
SU(2), = U ®U\Q)\1) , the optimal cloner is covariant and can be expressed 

as foUows, 

M 

|V',V'±) ^ ^a.-.MKM- j)^,^^) ® |(M- (5.34) 

3=0 

where the coefficients aj^M are given by 

\/3(M - 2j) 



Oij,M = {-ly 



+ 



^/2{M+í) ^/2M{M + 1){M + 2) 



(5.35) 



The cloning machine (5.34) is symmetric with respect to the interchange of \ip) 
and \ijj±). The cloner reqiüres an ancilla whose size is the same as the size of the 
output Hilbert space, i.e. the ancilla Hilbert space is also a symmetric subspace of 
M qubits. The ancilla contains M approximate copies of the state \tjj±) and the 
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Figuro 11: Cloning of a pair of orthogonal qubits by means of stimulated parametric 
down-conversion. BS - auxiliary beam splitters, Cj - nonlinear crystals, WP - 
wave plates for input-state preparation, M - mirrors, PD - triggering single-photon 
detectors for conditional preparation of single-photon states. 



fidelity of these anti-clones is the same as the fidelity of the clones. The single-clone 
fidelity can be calculated as weighted average of the coefíicients ce'j^M' 

^-^W = E^Í7^4m, (5.36) 
and after a simple àlgebra we arrive at 

The fidelity monotonically decreases with the increasing nuniber of clones M and in 
the limit M ^ oo we recover the fidelity (5.33) of the optimal state estimation from 
|'0)I'0-l)· Upon comparing the fidelity F^{M) to the fidelity of the optimal cloner for 
a pair of identical qubits, F\\{M) = (3Aí + 2)/(4M), we see that F\\{M) > Fj_{M) 
for M < 6, while F±{M) > -Fi|(M) for M > 6 and the cloner (5.34) outperforms 
the Standard universal cloner. 

We have seen in Section 5.1 that the optimal universal cloning of polariza- 
tion states of photons can be realized by means of stimulated parametric down- 
conversion. It turns out that also the optimal cloning with a pair of orthogonal 
qubits as the input can be performcd in the same way, if the photons in states \íp) 
and IV'-l) are fcd to the input signal and idler ports of the amplifier, respectively, 
as schematically illustrated in figure 11. We assume that the parametric amplifica- 
tion in the nonlinear crystal C3 is governed by the singlet-type Hamiltonian (5.1), 
which is invariant under the simultaneous rotation of the signal and idler qubits, 
(C/® U)H{W ® U^) = H. Assuming the input state |l)ay|0)aH|0)6y|l)b//, the 
output state after the amplification in the crystal C3 reads, 

00 

l^out) - ^""""'(1 - r^)!*^^,), (5.38) 

M=0 
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where the state with M clones and M anti dones is given by 

M 

|M/±,m) =^(-1)^' [(M-j)(1-A2)-A2] \M-3),v\j)aH\3)w\M-3)bH, (5.39) 

A = tanh(K;í) and í is the effective interaction time. 

In contrast to universal cloning with idcntical pliotons at tlic inpnt, the statc 
|Í'_l,m) depends on the strength of the parametric amphfication A. The cloner that 
produces M copies is obtained by post-selecting only the events when exactly M 
photons arc dctcctcd in signal and idlcr spatial modes, which corresponds to tlic 
selection of the state |\1/_l,m) i:roni the superposition (5.38). The fidelity of the 
cloner depends on A, 

FAM y) = 3.^-2.(2M+l) + |M(M+l) 

'^^ 6y2-6My + M(2M + l) ' ^^^^> 

where y = A^/(l — A^) = sinh^(Kí). The optimal parametric gain which maximizes 
the fidelity (5.40) can be found by solving the equation ^^-^^^'^^ = q which yields 



M 1 M{M + 2) 

On inserting the optimal y into Eq. (5.40) we recover the fidelity (5.37) hence 
the optimal cloning of a pair of orthogonal qubits can be achieved by means of 
stimulated parametric down-conversion with properly chosen gain. 



§ 6. Phase-covariant cloning of photons 

Up to now we have studied universal cloning machines that clone equally well all 
States. In many situations, however, one deals only with a subset of the states. An 
archctypal examplc is the Bennett-Brassard 1984 protocol for quantum key distri- 
bution (Bennett and Brassard, [1984]), which utilizes four non-orthogonal states 
|0), |1), |0) + |1), and |0) — |1). If we restrict the range of the admissible input 
states of the cloning machine, then wc can expcct that the machinc will exhibit 
better performance than the universal cloner and will reach higher fidelity. In this 
section we shall study the phase-covariant cloning machines which optimally clone 
all states that are balanced superpositions of the computational basis states, 

where the phases (pj can be arbitrary. 



6.1. PHASE-COVARIANT CLONING OF QUBITS 

The simplcist and perhaps most important example is the 1^2 phase-covariant 
cloning machine which can be used as the optimal individual eavesdropping attack 
on the BB84 protocol (Fuchs, Gisin, Griflaths, Niu and Peres [1997]). In contrast to 
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the universal cloners, the optimal cloning transformation depends on whether the 
single-clone fidelity or the global fidelity are taken as a figure of mèrit that should be 
maximized. In the context of eavesdropping on quantum key distribution protocol, 
it is natural to considcr the singlc-clonc fidclitics, since they quantify the amount 
of information transmittcd to the receivcr and gaincd by the eavcsdropper. 

The optimal symmetric 1^2 cloning transformation for qubits that maximizes 
the single-clonc fidelity has the foUowing form (Bruss, Cinchctti, D'Ariano and 
Macchiavello [2000], Bruss and Macchiavello [2001], Fan, Weihs, Matsumoto and 
Imai [2002]), 

|0)|Ai„) ^ -L|0)a|0)b|0)c + 1{\0)a\1)b + \1)a\0)b)\1)c, 

\1)\A^) - -^|1)^11)b|1)c + ^(|0)a|1)s + |1)a|0)b)|0)c. (6.2) 

The two clones are stored in qubits A and B and the fidelity of each clone reads 
Fi_^2 = (1 + l/v^)/2 ~ 0.855 which is indecd slightly higher than the fidelity 
of the optimal universal 1^2 cloner for qubits, F^^2 = 5/6 ~ 0.833. Note 
that besides a blank copy qubit, the transformation also requires another ancilla 
qubit C. Howcvcr. in contrast to universal cloning, this ancilla is not ncccssary 
and one can design a simplified cloning transformation which achieves the same 
fidelity and requires only two qubits: the input and a blank copy (Niu and GrifRths 
[1999], Durt and Du [2004]). This is very important from the experimental point of 
view since it is much easier to realize a two-qubit transformation than three qubit 
transformation. The econòmic phase covariant cloner can be obtained by projecting 
the ancilla C on the basis state ]0) (or ]1)). If we project on ]0), then we get 

|0)a|0)b - |0)a|0)b, |1)a|0)b-^(10)a|1)b + |1)a|0)b). (6.3) 

An alternative econòmic cloning transformation can be obtained from (6.3) by 
exchanging and 1. Interestingly, the cloning machine (6.3) is optimal not only 

for the States on the cquator of the Poincaré sphere but also for all the statcs on 
the northern hemisphere, i.e. aU states cos(é'/2)l0) + e'"^ sin(é'/2)ll) with 9 < Tr/2 
(Fiuràsek [2003]). The optimal asymmetric cloning machine which produces two 
clones with diffcrcnt fidclitics Fa and Fb is obtained by breaking the symmetry in 
the output superposition of )10) and )01), 

10)a10)b ^ 10)a10)b, ll)A|0)B^cos^lO)All)B + sint?]l)AlO)B, (6.4) 
and the two fidclitics can be expressed as follows. 



1 , „ 1 



Fa = -{í+sin^), Fb = -(1 + C0S1?), i? e 0,- . (6.5) 



2' " "2 



TT 



2 



The phase-covariant cloning machine that maximizes the global two-qubit fidelity 
has qualitatively similar structure as the cloner (6.2), 

10)1A.) - -L[10)^10)b|0)c + (|0)^11)b + |1)a|0)b)|1)c], 

\1)\A^) - ^[\1)a\1)b\í)c + mA\l)B + 11)a|0)b)|0)c], (6.6) 
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and it rcachcs fidclity Ff^2 =3/4 which is again higher than the global fidelity of 
the universal cloner, F^^^^''^ = 2/3. 

The phase-covariant cloning can be extended to the case when we posses N copies 
of the statc and would likc to preparo M cloncs, M > N. The optimal 1 ^ M 
phase-covariant cloning machine was determined by Fan, Matsumoto, Wang and 
Wadati [2001], who considered the single-clone fidelity as the figure of mèrit. The 
structure of the cloning transformation dcpcnds of the parity of M. If M is even, 
there exist two independent cloning transformations, 

|0)^|M,M/2-l), |l)^|M,M/2), (6.7) 

and 

|G) ^ |M, M/2), |1) ^ |M, M/2 + 1). (6.8) 

Note that there are in fact infinitely many cloning transformations since any convex 
mixture of the operations (6.7) and (6.8) is also optimal. On the other hand, if M 
is odd then we get only one optimal transformation 

|0)^ |M,(M-l)/2), |l)^|M,(M + l)/2). (6.9) 

The resulting fidelity is 



(6.10) 



The optimality of the cloning transformations (6.7), (6.8) and (6.9) can be proved 
using the same method that was employed in Sec. 4.2 to prové the optimality of 
the í ^ M universal cloning machine for qubits. In particular, the; singlc-clone 
fidelity can be expressed as F = Tr[SR], where S is the operator isomorphic to the 
cloning CP map and 

1 1 ^"^ 

ü= -70n+,M + ^(|O)(l| + |l)(O|)® ^ DM,k{\M,k + l){M,k\ + \M,k){M,k + l\), 

k=0 

(6.11) 

where DM.k = \/ (M — k){k + l)/M. The fidelity is uppcr bounded by the maxi- 
mum eigenvalue rmax oi R, F < 2rmax, and this bound is saturated by the above 
phase covar iant cloner s. 

Fan, Matsumoto, Wang and Wadati [2001] also conjecturcd the structure of the 
general optimal N M phase covariant cloning transformation for qubits. The 
proposed generalization is very straightforward, namely, every input symmetric A''- 
qubit statc \N, k) is transformcd to an M qubit symmetric state |M, k + j) with 
the constant j adjusted such that the fidelity is maximized. If A^ and M have the 
same parity, M = N + 2L, then the suggested cloning map is |A'',_j) — > \M,j + L), 
and the corresponding fidelity is 
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When the parity oi M and N difFers, M ~ A^ + 2L + 1, thcn the two possible cloning 
transformations are either |iV, j) |AÍ, j + L), or \N,j) \M,j + L + 1), and the 
corresponding fidelity is 

F'n-^m = 2 + M^g\/(^)G + l) 

X [ V(iV + L~j + l){L + j + l) + ^iL + j+2){N + L- j)]. 

(6.13) 

The optimality of the fidelity (6.12) was proved in D'Ariano and Macchiavello [2003] 
by exploiting the genèric theory of covariant cloning machines, see D'Ariano and Lo 
Presti [2001]. In contrast, if N and M have different parities, the optimal phasc- 
covariant cloning transformation that was found in D'Ariano and Macchiavello 
[2003] difïers from Eq. (6.13). 



6.2. PHASE-COVARIANT CLONING OF QUDITS 

Going beyond the cloning of qubits, the 1^2 phase covariant cloning of qudits 
(6.1) has been investigated (Fan, Imai, Matsumoto and Wang [2003], Lamoureux 
and Cerf [2005], Rczakhani, Siadatncjad and Ghadcri [2005]). It can be shown that 
the optimal cloning transformation for qudits (6.1) has the structure 

d-i 

\j) - a\jj)AB\j)c + -7==^, Y.^\jI)ab + \13)ab)\1)c, (6.14) 
V2(CÍ - 1) i^- 

where 0^+0^ = 1. The two clones are contained in qudits A and B while the qudit 

C serves as an ancilla. Notc that Eq. (6.14) is a direct extcnsion of the cloning 
transformation for qubits (6.3). The coefScients a and /3 have to be optimized such 
that the cloning fidelity is maximized. After some àlgebra one arrives at 

/l d-2 ^ (\ d-2 

"=l2- 2Vd^+4c^-J ' ^=U + 2Vd^+4d-J ' ^'-''^ 

and the fidelity reads 

In contrast to the phase covariant cloning of qubits, we cannot get rid of the an- 
cilla C because if we project ancilla on the computational basis state \k) then the 
conditional map is not unitary. So for d > 2 it seems impossible to implement the 
optimal phase covariant 1 — » 2 cloning in an econòmic way, without ancilla. 



6.3. OPTICAL PHASE-COVARIANT CLONING 

In contrast to the universal cloning, the optical experimental implementation of 
phase covariant cloning machines has received much less attention. This may come 
as a surprise in view of the apparent simplicity of the optimal cloning transforma- 
tion (6.3). However, the phase-covariant cloning exhibits much less symmetry than 
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Figure 12: Phase-covariant cloning using an interference of two phhotons on an un- 
balanced beam splitter. (a) Scheme with a single beam splitter BS that difFerently 
reflects vertical and horizontal polarizations. (b) An alternative setup involving po- 
larizing beam splitter PBS and an unbalanced beam splitter BS' whose reflectance 
does not depend on the polarization. 

the universal copying and the methods such as stimulated amplification or sym- 
metrization cannot be readily extended to implement the 1 — > 2 phase covariant 
cloning machine. 

It is nevertheless possible to conditionally realize the 1-^2 phase covariant 
cloning of photonic qubits with linear òptics (Fiuràsek [2003] ) . As usual, the qubits 
are encoded into polarization states of single photons, and the state to bc cloned 
is a balanced superposition of vertical and horizontal polarization, \ip) — {\V) + 
é^'^\H)) / ^/2. Besides the input state, the cloning requires also a second photon, the 
blank copy which we assume to be initially prepared in the state \V). Written in 
the basis of polarization states, the cloning transformation (6.3) becomes 



The scheme of the cloning machine is illustrated in Fig. 12(a). The input photon 
and the blank copy are combined on an unbalanced beam splitter whose trans- 
mittance tj and reflectance Vj for the vertical [j = V) and horizontal [j — H) 
polarizations are different. Only the events when the two photons leave the beam 
splitter in different output ports are post-selected. The principle of operation of 
the cloner is easy to grasp. If the input is in state \V), the two photons at the 
output must be in state \VV) since the blank copy is initially in the state \V). On 
the other hand, if the input to be cloned would be in the state \H) then the beam 
splitter would produce a superposition of \HV) and \ VH). By properly choosing Vj 
this superposition can be made balanced and the conditional map becomes exactly 
the unitary (6.17). 

The mixing of the modes on a beam splitter is governed by the linear canonical 
tr ansformations : 



only the events when there is one photon in the left output arm (mode a) and one 
photon in the right output arm (mode b) reads 



\VV)^{rl-tl)\VV), \HV) ^rHrv\HV)-tHtv\VH). (6.19) 



\VV) ^ \VV), 



\HV)^^{\HV) + \VH)). 



(6.17) 
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This transformation becomes fuUy equivalent to Eq. (6.17) if the following condi- 
tions are satisfied, 

4 - = ^rarv = -V^íhív (6.20) 

These constraints iniply that rn = ty, tn = — ry and (ry — ty) = \f2r\'ty . On 
combining this equation with the normahzation ry + íy = 1 we can determine ry. 
After simple àlgebra we obtain 



2 V \/3 

The probability of successful realization of the phase-covariant cloner is given by 

P={T\-t\f = \. (6.22) 

The required beam splitter with different transmittances for vertical and hori- 

zontal polarizations can bc simulatcd by a Mach-Zchndcr intcrfcromctcr with the 
polarization dependent phase shifters in its arms, such as Soleil-Babinet compen- 
sators, so that the phase shift and, consequently, the splitting ratio could be con- 
troUcd indcpcndcntly for vertical and horizontal polarizations. The sctup could bc 
also modified to work with a beam splitter whose reflectance is the same for both 
vertical and horizontal polarizations. This alternative configuration is depicted in 
Fig. 12(b). The signal and blank copy photons arc first combincd on a polarizing 
beam splitter PBS that reflects vertically polarized photons and transmits horizon- 
tally polarized photons. The two beams are then re-combined on a beam splitter 
with reflectance r. If the signal photon is initially vertically polarized, thcn a verti- 
cally polarized photon enters each input port of BS. If the signal photon is polarized 
horizontally, then it is switched to the right arm and two photons in orthogonal 
polarization statcs impingc on the right input port of BS. The polarizing beam split- 
ter ensures that the role of the transmittance and reflectance for the horizontally 
polarized photon is interchanged with respect to the scheme shown in Fig. 12(a). 
It can bc easily shown that this sctup Icads to the cloning transformation (6.17) 
provided that the reflectance of the beam splitter is equal to r-^ = (1 -|- l/·\/3)/2. 

In the experiment, it may not be easy to precisely control the transmittance. It is 
thcrcforc important to invcstigatc how the pcrformancc of the sctup shown in Fig. 
12(b) depends on the reflectance r of the beam splitter. The cloning transformation 
remains phase-covariant and the cloning fidelity F is the same for both clones and 
docs not dcpcnd on 0. Howcver, F bccomcs a function of r. After some àlgebra 
one arrives at the formula for the fidelity of cloning of equatorial qubits, 



F=^- 
2 



^ 2r(2r^ - 1)7(1^ 
2r4 - 2r2 + 1 



(6.23) 



It turns out that the cloning is rather robust with respect to the variations of the 
reflectance of the beam splitter and cloning fidelity F > 0.8 can be achieved for a 
broad range of beam splitter reflectances 0.7 < r^ < 0.9. 

In the experiment, the required pair of photons can be produced in spontaneous 
Type-I parametric down-conversion and the dcsircd initial states of the photons 
can be prepared with the usc of wavc plates. After cloning, the states of the two 
clones can be analyzed by a sequence of wave plates, polarizing beam splitters, and 
single photon detectors, similar ly as in the experiments on universal cloning. 
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6.4. EXPERIMENTAL l-TO-3 PH ASE- CO VARIANT CLONING 

Remarkably, while the optimal 1—^2 phase covariant cloning transformation (6.2) 
has not yct bccn implemented for optical qubits, thc optimal 1^3 phasc-covariant 
cloning of polarization state of a single photon has been demonstrated experimen- 
tally by Sciarrino and De Martini [2004]. The set of cloned states included all 
linear polarization states cos9\V) + siní^jiï). Thc first stcp in thc copying proccss 
consisted of the optimal 1 — » 2 universal cloner describcd in Section 5.1 which pro- 
duced two clones and one anti-clone. In the next stcp, the anti-clone was converted 
into a clonc by applying a unitary transformation o"y with thc; hclp of a lialf-wavc 
plate. The final step was to symmetrize the state of the three clones by combining 
the two cones and the anti-clone on a balanced beam splitter and selecting only the 
evcnts whcrc all thrcc photons cndcd up in thc same output spatial mode. In this 
way thrcc copies of cqual fidelity wcrc produccd. 

The experimentally observed fidelities were -Fí!íl3(|+)) = 0.76 for the state |+) = 
2-V2(|y) + \H)) and Fl^^{\H)) = 0.80. This should be compared with the theo- 
rctical maximum Ff^r^ = 5/6 ~ 0.833. It is also instructivc to makc a comparison 
with the fidelity of the optimal universal 1^3 cloner, F^I^'^ = 7/9 0.778. One 
can conclude that the experimental phase- covariant cloning machine operates very 
close to its thcorctical limit and for ccrtain inputs it achieves better fidelity than 
what would be possible with universal cloning machine. 



§ 7. Cloning of optical continuous variables 

So far we have considered cloning of quantum states in finite dimensional Hilbert 
spaccs. Diiring recent ycars, howcvcr, quantum Information proccssing in systems 
with infinite dimensional Hilbert space, such as modes of the electromagnètic field, 
has attracted great deal of attention, see, e.g, Braunstein and van Loock [2005] . In 
this approach, thc quantum information is usually cncodcd into two noncommut- 
ing quadrature operators x and p which satisfy canonical commutation relations 
[x,p] = i. Since these operators have continuous spectra, one speaks about quan- 
tum information processing with continuous variables. 

The universal cloning machine for states belonging to infinite dimensional Hilbert 
space can be formally obtained as a limit of the universal cloning machine for qudits 
when d ^ 00. One finds that the single-clone fidelity of the universal 1^2 cloner 
is 1/2 which means that the optimal cloning can be achicved be a very simple 
strategy where the input state is sent with probability 1/2 to the first or second 
output, while the other output is prepared in maximally mixed state. Besides being 
rather trivial, this universal cloner is not of great practical interest because most 
of the quantum information protocols with continuous variables involve only the 
so-called Gaussian states. These state have Gaussian Wigner function and their 
great advantagc is that thcy can be rclativcly easily generatcd and manipulated 
in the lab with the help of linear optical interferometers and optical parametric 
amplifiers which produce squeezed and entangled Gaussian states. 
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7.1. CLONING OF COHERENT STATES 

Among the Gaussian states, the coherent state is perhaps the best known exam- 
ple. The coherent state |a) can be defined as a displaced vacuum state D(q;)|0), 
where D{a) = exp{aa^ — Q*a) is the displacement operator. Coherent state is 
the eigcnstatc of the annihilation operator, a\a) = a\a) and it is also a minimum 
uncertainty state. The variance of all rotated quadratures xe = x cos + p sin 9 is 
the same and equal to 1/2. The Glauber P-distribution of the coherent state is a 
Dirac delta function. In this sense, the coherent states are not usually considered 
as non-classical states in the quantum-optical sense. Still, these are purc quantum 
states and they carry quantum noise. This makes these states suitable for applica- 
tions siich as quantum key distribution. It has bcen shown thcorctically (Grosshans 
and Grangier [2002], Grosshans, Cerf, Wenger, Tualle-Brouri and Grangicr [2003], 
Grosshans and Ccrf [2004], Iblisdir, Van Assche and Cerf [2004]) and dcmonstrated 
experimcntally (Grosshans, Van Assche, Wenger, Tualle-Brouri, Ccrf and Grang- 
ier [2003]) that secure key distribution can be achieved with coherent states and 
balanced homodyne detcction. 

Let us first consider the optimal Gaussian cloning of coherent states introduced 
in Cerf, Ipc and Rottenberg [2000], Lindblad [2000]. In this scenario, the class 
of admissible cloning transformations is restricted to Gaussian operations, which 
preservo the Gaussian shape of the Wigner function. Intuitivcly, one could cxpect 
that the Gaussian cloning should be optimal. This is indced truc if the figuro of 
mèrit is the global M-clone fidelity or if the quality of the clonos is quantifiod in 
terms of the noise added to the two quadratures x and p (Cerf and Iblisdir [2000]). 
However, it has becn realized reccntly that, remarkably, the single-clone fidelity of 
the 1 — > 2 cloning of coherent states is maximized by a non-Gaussian cloner (Cerf, 
Krüeger, Navez, Werner and Wolf [2005]). 

Let us begin with the Gaussian 1^2 cloning. We require that the mean vàlues 
of the quadratures of the two clones A and B are equal to the mean vàlues of 
the quadratures x-m and pin of the input coherent state a. This guarantees that 
the cloning transformation is invariant with respect to the displacements and the 
cloning fidelity does not depend on the amplitude a. In the Heisenberg picture, the 
most general Gaussian cloning transformation can be written in the form, 

XA = Xí^-\- XA, Xb = Xí^-\- Xb, (7.1) 
PA=Pi-a+PA, PB =Pin+PB- (7.2) 

The operators xa, Pa, xb, Pb represent the noise that is added to the two copies 
during the cloning process and they all commute with a;i„ and Pi^. The quadrature 
operators must satisfy the canonical commutation relations, which implies 

[xa,Pb] = -i, [xb,Pa] = -i- (7.3) 

The Heisenberg uncertainty relation gives a lower bound on the products of the 
variances of the noise operators, 

{{AxAmApBf) > l, {{AxBmApAf) > \. (7.4) 

The cloning should add noise isotropically, the variance of the x and p quadratures 
of each clone should be the same. Since the noise operators are not correlated with 
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Xin and pin, thc variance of thc quadratures (7.2) of the two clones is sum of two 
variances, and the isotropy condition is satisficd if 

{{AxAf) = {{ApAr) = ÜA, ((AÍb)') = {{ApB?) = UB. (7.5) 

The two uncertainty relations (7.4) the boil down to a single constraint 

ÜAnB > ^- (7.6) 

The state of each clone is a mixed Gaussian state, namely a coherent state with 
added thermal noise with mean number of thermal photons equal to üj, j = A, B. 
The fidehty of cloning can be most easily calculated from the Husimi Q-function, 
which is defined as the overlap of the density matrix with coherent state, Q(/3) = 
7r~^(/3|p|/3). The Q function of the j-th clone reads (Fiuràsek [2001a]), 



1 + n,- 



(7.7) 



Thc fidehty can be calculated as Fj{a) = TTQj{a) = 1/ {1 + üj). Thc best trade-off 
between the fidelities of the two clones is obtained when the equality holds in Eq. 
(7.6), and we get 

= ^> = (7.8) 

where 7 is a parameter which controls the asymmetry of the cloning. The fidelity 
of the optimal 1 — > 2 symmetric (7 = 0) Gaussian cloner is F = 2/3. 

7.2. CLONING BY PHASE-INSENSITIVE AMPLIFICATION 

If thc coherent states arc carried by optical modes, then the cloning can be real- 
ized with the usc of a phasc insensitive amplification of light (Braunstein, Cerf, 
Iblisdir, van Loock and Massar [2001], Fiuràsek [2001a]). This is very natural and 
intuitive result, because the idealized perfect cloning amounts to noiseless ampli- 
fication of the coherent state, \a) \^/2a). The optimal amplification that adds 
the minimum amount of noise can be performed, e.g., in nondegenerate optical 
parametric amplifier (NOPA), which transforms the input annihilation operator a 

as ttout 

= VGüin + VG — Icjjj, where c is the annihilation operator of the idler 
mode in NOPA. The setup for asymmetric cloning of coherent states is shown in 
Fig. 13(a). It consists of a Mach-Zehnder interferometer with an amplifier in one of 
its arms. The signal is initially divided into two beams and one beam is amplified 
such that the total mcan intensity is twice the input intensity. The two clones are 
obtained by re-combining the two beams on the second beam splitter. The splitting 
ratios of the unbalanced beam splitters BSi and BS2 and the intensity gain G of 
the amplifier can be expressed in terms of the asymmetry parameter 7 as follows, 



\/2 sinh 7 



27 



G = 1 + cosh(27), Í2 = . (7.9) 



Vl + 2sinh''7 VI + 

The setup becomes particularly simple for symmetric cloner. In this case the first 
beam splitter disappears and the whole input signal is amplified with gain G = 2 
and then divided into two modes on a balanced beam splitter, see Fig. 13(b). 
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Figure 13: Cloning of coherent states in a non-degenerate optical parametric ampli- 
íier. (a) Asymmetric 1 — > 2 cloner which consists of a Mach-Zhender interferometer 
with a non-degenerate parametric amplifier (NOPA) in one of its arms. The am- 
pHfication gain and the sphtting ratios determine the asymmetry of the cloner. (b) 
Simphfied scheme of a symmetric cloner. 



The procedure for symmetric cloning can be readily extended to the optimal 
symmetric N ^ M Gaussian cloning of coherent states. The cloning consists of 
three steps. First, the whole signal is coUected in a single mode using an array 
of — 1 beam splitters with properly chosen transmittances, \a)®^ \^/Na). 
Next, the collected signal is amplified with a gain G = M/N. Finally, the amplified 
signal is distributed among the M modes with the help of another array of M — 1 
unbalanced beam splitters such that the mean complex amplitude in each mode is 
the same and equal to a. The fidelity of this cloner does not depend on a and each 
clone is in a coherent state with thermal noise described by the Husimi function 
(7.7). The total mean number of thermal photons in all modes is G — 1 = M/N — 1 
and the noise is equally divided into M modes hence the thermal noise in each 
clone is fi = 1/A^ — 1/M. On inserting this into the expression for the fidelity, 
F = 1/(1 + n), we obtain 

. (7.10) 

MN + M -N ^ ' 

In the limit of infinite number of copies, M — > oo, we get F — N / [N + 1) which is 
the fidelity of optimal estimation of a coherent state from A^ copies. Similarly as in 
the case of universal cloning of qubits (Bruss, Ekert and Macchiavello [1998]), the 
connection between optimal cloning and optimal state estimation can be exploited 
to prové that (7.10) is the maximal fidelity of the Gaussian A^ — ^ M cloning of 
coherent states (Cerf and Iblisdir [2000]). As shown by van Loock and Braunstein 
[2001], it is also possible to clone coherent states via an extended continuous- variable 
teleportation. The telecloning requires a specific multimode entangled Gaussian 
state that can be generated by mixing single-mode squeezed vacuum states on an 
array of unbalanced beam splitters (van Loock and Braunstein [2000]). 

As already noted before, the Gaussian machine depicted in Fig. 13 is not the 
optimal one if the single-clonc fidelity is taken as a figure of mèrit (Cerf, Krüeger, 
Navez, Werncr and Wolf [2005]). The optimal non-Gaussian cloner can achieve a fi- 
delity Fijiax = 0.6826, which is slightly higher than the maximum fidelity achievable 
by Gaussian transformations, F = 2/3 « 0.6667. Interestingly, the optimal non- 
Gaussian cloner can be obtained from the setup shown in Fig. 13(b) if the input 
ports of the idler mode of the amplifier and the auxiliary mode of the beam splitter 
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Figure 14: Cloning of coherent states using only linear òptics and homodyne de- 
tection. í^i, U2 and vy, denote vacuum auxiliary input modes, D is a displacement 
proportional to the measurement A. (After Andersen, Josse and Leuchs [2005].) 



are fed with a specific non-Gaussian entangled state \iI})bc = Y1^=q c„|2n)s|2n)c- 
The coefíicients Cn can be optimized in order to maximize the cloning fidelity which 
results in the above given value F^ax = 0.6826. It should be stressed that while 
the non-Gaussian cloner maximizes the fidelity, the variance of the quadratures of 
the clones is higher than for the optimal Gaussian cloner. In applications such as 
quantum key distribution with coherent states and balanced homodyning, where 
the quantum channcl bctween Alice and Bob is characterized in terms of the first 
and second moments of the transmitted quadratures, the aim of the eavesdropper 
is to minimize the quadrature variance instead of the fidelity and the Gaussian 
cloning (or its variant entangling Gaussian cloner in case of reverse reconciliation 
protocol) can be the most dangerous individual eavesdropping attack. 

7.3. EXPERIMENTAL CLONING OF COHERENT STATES 

Recently, optimal Gaussian 1 2 cloning of coherent states has been experimen- 
tally demonstrated by Andersen, Josse and Leuchs [2005]. The distinct feature of 
this experiment is that it docs not require an amplifier which is replaced by a clever 
combination of measurement and feedback. A simplified scheme of the experimental 
setup is depicted in Fig. 14. The mode A contains the coherent state to be cloned. 
The beam is split into two parts on a balanced beam splitter whose auxiliary input 
port vi is in vacuum state. The output annihilation operators thus read 

a' = + í^l.in), v'l = ~ í^l.in)- 

The output beam v'^ is sent to an eight-port homodyne detector, which consists 
of a balanced beam splitter foUowed by two balanced homodyne detectors. This 
detector efïectively measures the operator A = v[ + v\, where v\ is the creation 
operator of an auxiliary vacuum mode. After the measurement, the mode a' is 
displaced by the amount A which is in practice achieved by mixing this beam with 
a strong coherent beam with amplitude X/y/T--T on a highly unbalanced beam 
splitter with transmittancc T « 99 %. The resulting displaced beam is effectively 
the amplified input, 

Odisp = \/2ain + i^l ■ 
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The cloning is finished by dividing the amplificd bcam into two parts with the help 
of another balanced beam splitter, thereby preparing the two clones of the input 
coherent state. The fidehty observed in the experiment was about 65% which is 
very close to the optimal value 2/3 ~ 0.667. 

7.4. GAUSSIAN DISTRIBUTION WITH FINITE WIDTH 

Up to now, we have assumed that the distribution of the coherent states that 
should be cloned is uniform over the whole phasc spacc. Howcvcr, this is clearly 
an idealization, since the mean energy of the input state would be infinite. A more 
reaHstic scenario, considered by Cochrane, Ralph and Dohnska [2004], is that the 
coherent states are drawn from a Gaussian distribution with width a and centered 
on vacuum, so that the a-priori probabiüty that the cloned state is \a) is given by 

1 / 



This occurs for instancc in the quantum kcy distribution with coherent states 
(Grosshans, Van Assche, Wenger, TuaUe-Brouri, Cerf and Grangier [2003]). 

If the width a of the distribution (7.11) is finite, then we possess some information 
that can be explorcd in ordcr to incrcasc the avcragc cloning fidcHty. Also. the 
probability (7.11) is not invariant with respect to the displacements, so there is no 
reason to search for covariant cloner. The fidelity of the cloner may depend on the 
input state and the figure of mèrit that should be maximized is the average fidelity, 

T = j P{a)F{a)d?a. (7.12) 

It turns out that the optimal symnictric; 1^2 Gaussian cloning transformation is 
still the amplification foUowed by the beam splitting on a balanced beam splitter. 
However, the gain G depends on a. After the amplification and beam splitting, 
the coherent amplitudc in cach mode is ayJG/2 and the mcan numbcr of chaotic 
photons in each mode is rï = (G — 1)/2. The fidelity of cloning a particular coherent 
state \a) reads 

2(1-7^/2)2 



G + l 



(7.13) 



After the averaging over the Gaussian distribution (7.11) we arrive at the expression 
for the mean fidelity, 

T= ? =. (7.14) 

G+ 1 + 2(72(2 + G-2V2G) 

We have to find the maximum of T undcr the constraint G > 1. It turns out that 
there are two different solutions in dependence on the value of cr. If a'^ > a^j^ = 
(1 + -\/2)/2, then it is optimal to amplify the signal and the optimal gain is 
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On the othcr hand, if cr^ < cr^j^ thcn it is optinial to simply dividc thc input signal 
into two beams without any aniplification, and G — 1. The resulting cloning fidelity 
is 

4cr2 + 2 



6cr2 + l' 
1 



2 \ 2 



(7.16) 



a < a. 



th- 



1 + (3 - 2\/2)ct2 ' 

The average fidehty monotonically increases with the decreasing width of the dis- 
tribution (7.11) and in the limit a ^ we get !F = 1, as expected. 



7.5. CLONING OF CONJUGATE COHERENT STATES 

In Scction 5.4 wc have discussed a cloning machinc for a pair of orthogonal qnbits. 
This device possesses a natural and very interesting continuous- variable analogue, 
namely, one can consider a cloning machine for coherent states \a) whose input 
consists of N copies of thc statc \a) and of N' copies of thc complex conjugatc 
coherent state \a*). This problem was analyzed in detall in even more general 
setting in Cerf and Iblisdir [2001a]. 

Without any loss of generality, wc can assunic that a pair of arrays of bcam 
splitters is used to coUect all signal into two modes and the input state of the 
cloning machine thus reads \VNa)A\VN'a*)B- The goal of cloning is to produce 
M copies of \a) with minimum addcd noisc. This could be again accomplishcd with 
the help of the non-degenerate parametric amplifier. While the mode A represents 
the signal input similarly as before, the mode B is sent to the idler input port of the 
amplifier. Assuming amplification with intensity gain G, the output annihilation 
operator of the signal mode is given by, Oout = VGuin + — 1 b}^. Note that both 
terms Oin and in the above formula contribute to the total coherent signal in 
aout- If the cloning should be performed with unity gain, then G must satisfy 

Vm = VgVn + VG^VN', (7.17) 

and G can be easily determined by solving the above quadratic equation. 

A careful analysis reveals that for certain vàlues of N, N' and M the cloning 
with conjugate inputs could be more efficient than the Standard cloning of coherent 
states. To be fair, we should compare the cloning fidelities for inputs consisting 
either of iV + A^' copies of |a) or including N copies of |q;) and N' copies of |a*). 
The advantage of dealing with complex conjugate inputs could be most easily il- 
lustrated in the limit of infinitc number of clones, M oo, wherc the optimal 
cloning becomcs equivalent with optimal state estimation. It has becn shown in 
Cerf and Iblisdir [2001b] that when possessing a single copy of ja) ja*) we can es- 
timate \a) with fidelity Fc.c. = 4/5 which is strictly higher than the estimation 
fidelity F — 2/3 corresponding to the input state la)'^^. In the former case the 
optimal dctcction strategy is the nonlocal continuous-variable Bell measurement 
where the quadratures xa+ xb and pa— Pb are measured simultaneously. 



§ 8. Conclusions 



The quantum no-cloning theorem is a crucial aspect of modern quantum mechanics 
and one of the cornerstones of quantum Information theory. Besides its fundamen- 
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tal intcrcst for thc foimdations of quantum physics, the impossibility of exactly 
copying an unknown quantum state is crucial for the security of the quantum key 
distribution protocols. Going beyond the no-cloning theorem, it is possible to de- 
sign approximate quantum cloning machines, which enable the copying of quantum 
Information in an optimal - albeit imperfect - way, an issue which has attracted a 
considerable attention over the last decade. 

This review paper aims at providing an exhaustive view of the various quantum 
cloning machines that have becn introduced since this concept was put forward 
by Buzck and Hillery in 1996. The mathematical description of quantum cloning 
machines bascd on the isomorphism bctwecn maps and operators is dcvclopcd in 
detalls. A special attention is also paid to the experimental optical implcmcntations 
of these machines. The cloning of single photons has been accomplished by several 
groups by now, and these experiments represent a very valuable contribution to the 
toolbox of available optical methods for quantum Information processing. 

In the coursc of ycars. quantum cloning has grown into a genuine subfield of quan- 
tum Information scicnces, which is still currently very active both on the theoretical 
and experimental sides. The advanced methods of preparation, manipulation, and 
measurement of quantum states of light, whose development has been stimulated 
to a large extent by the perspectives of quantum Information processing, have re- 
cently enabled the demonstration of even more complex cloning machines. In the 
years to come, we anticipate many new achievements and breakthroughs in quan- 
tum Information sciences, and there is no doubt that quantum cloning will play an 
important role in these future developments. 
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